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ABSTRACT 


The mass, momentum, and energy flows through a quasi- 
steady, self-field MPD accelerator are measured over a 
current range of 8 to 50 kiloamperes and inlet mass flows 
of 2 to 36 grams per second of argon. Time-resolved arc 
chamber pressure and impact pressure in the exhaust jet, 
measured with a piezoelectric pressure probe capable of 
operation with a flat response from 1 to 140 kilocycles, 
are used to calculate the momentum flux profiles, the gas- 
dynamic contribution to thrust, and the total thrust. The 
momentum flux profile' indicates that the accelerator pro- 
duces a uniform, 2-inch diameter axial jet at the anode 
which expands into a Gaussian profile at an axial station 
11 inches from the anode. The electromagnetic component 
of the thrust is found to follow the familiar quadratic 
dependence on arc current, while a more complex empirical 
relation is needed to correlate the gasdynamic contribu- 
tion with the current and mass flow rate. Using available 
time-of-f light velocity profiles at a current of 16 kilo- 
amperes and a mass flow of 5.9 grams per second, calculated 
flux profiles of mass and kinetic energy exhibit a tendency 
for some fraction of the inlet mass flow to leak out at a 
low velocity around the central high velocity core. Dis- 
charge chamber geometry changes to eliminate the leakage 
will doubtless improve the arc performance. 
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CHAPTER 1 
INTRODUCTION 

Experimental investigations of the operation of plasma 
accelerators can be divided into two classifications. In 
one, the accelerator is treated as a black box, and the out- 
put thrust is measured as a function of the input mass flow 
rate and the electrical power. This information is then used 
to calculate the thruster performance (Ref. 1-1) . In the 
second type of investigation, the detailed structure of the 
electrical fields, magnetic fields and the thermodynamic state 
of the plasma are measured as a function of the input mass 
and energy and also as a function of time and position in the 
acceleration region. This more detailed information is used 
to study the acceleration processes and losses in the plasma 
accelerator (Ref. 1-2) . Both types of investigation have a 
common goal, the development of improved thrusters, but until 
recently, neither has been of much help to the other in achiev- 
ing that goal. 

The reason for this division has been the rapid evolution 
of plasma accelerators and their theory and measurement tech- 
niques. Each individual laboratory investigated its own 
unique operating conditions and accelerator geometry (eg. "T" 
gun, Faraday accelerator or Hall thruster) (Ref. 1-3, 4, 5, 6). 
Theoretical understanding of accelerator operation was insuf- 
ficient to allow reasonable comparisons of the widely disparate 
results, and frequently the measurement techniques were later 
found to be incorrect. Progress in finding better thruster 
geometries was rapid, and those who concentrated on improved 
detailed diagnostics and theoretical models were furnished 
with new geometries faster than they could be analyzed. At 
the same time, those who measured performance, and then went on 
to a new geometry found the cut and try method was much more 
successful and faster than relying on suggestions from the 
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analytically inclined laboratories. This led naturally to 
a deep division between the two types of experimental inves- 
tigations of plasma accelerators. 

The pace of change has slowed considerably since the 
invention of the MPD accelerator in 1964 (Ref. 1-7) . Most 
of the work on alternate geometries has been dropped. One 
of the three forms of the MPD accelerator (steady-state 
applied field, steady-state self field and pulsed self field, 
(Ref. 1-8) are studied in most plasma accelerator laboratories. 
No other radically different geometry or set of operating con- 
ditions has been found that provides a major improvement in 
performance. The emphasis has shifted to relatively minor 
geometry changes, such as the hollow cathode (Ref. 1-9), or 
to the search for the optimum value of operating parameters 
such as applied field strength or current-to-mass flow rate 
ratio (Ref. 1-10) . This has also allowed the more analytical- 
ly inclined investigators to make an increased contribution 
to plasma accelerator development, since further improvements 
must be based firmly on validated experimental measurements. 

The division between the two types of experimental investiga- 
tions is becoming less pronounced. 

The investigation described in the following chapters is 
a continuation of this process. Clark (Ref. 1-11) , Turchi 
(Ref. 1-12) and Oberth (Ref. 1-13) have measured the detailed 
structure of the electric and magnetic fields in the accelera- 
tion region of the Princeton MPD accelerator and used this 
information to study acceleration processes and losses. On 
this same accelerator, Clark (Ref. 1-11) has also measured 
input mass flow rates and electrical power to the accelerator 
and in addition, time-of-f light velocity profiles were ob- 
tained (Ref. 1-14) in the exhaust plume. The present investi- 
gation was designed to augment this information by providing 
a detailed picture of the mass, momentum and kinetic energy 
flow to and from the accelerator. This information, coupled 
with the detailed picture of the acceleration region, brings 
us one step nearer the goal of the analytical type investigation 
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gaining sufficient understanding of MPD acceleration so that 
improved thrusters can be built. At the same time, the new 
knowledge of the momentum flux can be used to calculate thrust. 
This completes the set of measurements required and used to 
evaluate performance: that is, thrust, inlet mass flow rate 
and electrical power,, The investigations described in the 
following chapters, therefore, provide a link between the 
performance and analytical types of investigation. 

The specific approach for determining mass, momentum 
and energy flow was determined from a control volume analysis. 
For the steady state time period (Ref. 1—11) , a control volume 
enclosing the entire acceleration region contains a constant 
amount of mass, momentum and energy. During this time, there- 
fore, the inlet and exit mass and energy flow rates must be 
equal and the surface and volume forces on the plasma must 
equal the rate at which momentum leaves the control volume. 

The dominant contributions to mass, momentum and energy input 
to and output from the control volume were estimated, and 
the measurements to be taken were selected to complete the 
picture. 

Figure 1-1 shows the control volume selected and the best 
estimate as to the dominant contributions to inlet and outlet 
mass, momentum and energy. The control volume was selected 
to make a mass, momentum and energy balance as easy as pos- 
sible. Its surface coincides with the electrode and in- 
sulator surfaces in the acceleration chamber and extends 
downstream sufficiently far so that all plasma currents are 
included and far enough to the sides to include the expanding 
plasma plume. With this control volume, there are no surface 
forces or flux through the sides of the volume and the entire 
j x B body forces and I V electrical energy is contained in 
the volume. 

The mass balance over this control volume was assumed 
to be dominated by argon mass flow supplied to the cham- 
ber and the mass flux in the accelerated jet. Ablated mass 
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was assumed to be negligible. Clark (Ref. 1-11) obtained 
approximate values for the inlet argon flow (hereafter re- 
ferred to as "nominal" mass flow rate) so only the mass flow 
in the jet was required to complete the mass balance. This 
mass balance was expected to show any errors in the "nominal" 
inlet flow or ablative mass (but did not because of mass 
leakage) . 

The principal axial force accelerating the plasma was 

.a jk 

assumed to be the j x B body force producing a net axial momen- 
tum flux in the exhaust plume. Since the field studies (Ref. 

-A 

1-11) showed axial symmetry, the j x B contribution could be 
calculated (Ref. 1-8) . Surface forces were assumed to be 
small except possibly for the axial pressure force exerted by 
the rear chamber wall on the plasma. In arc jets, the evolu- 
tionary ancestor of the MPD accelerator, this surface force 
was much greater than the integrated j x B body force and 
there is still some debate as to the relative contribution 
of these two forces. Therefore, to achieve a momentum balance 
over the control volume (and incidently thrust) , it is neces- 
sary to measure the momentum flow out of the downstream face 
of the control volume and also the upstream surface force 
(or gasdynamic contribution) . 

The principal energy flow into and out of the control 
volume were assumed to be the input electrical power, I V 
(Ref. 1-11), the cathode and anode losses (Ref. 1-12,13), and 
the kinetic and internal energy flows out of the downstream 
face of the control volume. The first three have been measured 
so if the kinetic and thermodynamic flux can be determined, an 
energy balance can be made. 

To complete a mass, momentum and energy balance, there- 
fore, measurements of the exit flux of these quantities plus the 
electrothermal thrust are required. Measurements of other dom- 
inant terms were already available. The outflow quantities to 
be measured are listed at the bottom of Fig. 1-1 and it is ap- 
parent that measurement of two properties, velocity and density 
would be sufficient to determine most of the outflow quantities. 
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Since velocity profiles were already available (Ref. 1-14) , 
only three additional measurements were required. These 
were: 1) electrothermal thrust; 2) frozen flow power, and 
3) some quantity related to the density at the downstream 
face of the control volume. It was decided to neglect frozen 
flow power and concentrate on developing and using a 
diagnostic tool for measuring the other two quantities. 

The diagnostic tool selected for development was a 
pressure probe. This development of a probe satisfying the 
environmental, sensitivity, time response and accuracy re- 
quirements, represents one of the major contributions of 
this investigation. This development is described in Chap- 
ter 2. The probe was used to provide a measurement of the 
electrothermal contribution of thrust by determining the 
plasma pressure in the acceleration chamber. The surface 
force on the plasma was then determined by multiplying this 
pressure times the surface area over which it acts. This 
use of the pressure probe and the results obtained are de- 
scribed in Chapter 4, The same probe could also be used to 
determine the flux of momentum, mass and kinetic energy from 
the downstream face of the control volume. Used as the pres- 
sure sensor in an impact (pitot) probe, the quantity yOu 
could be measured as a function of radial position. By in- 
tegrating over the downstream face of the control volume, 
the momentum flow (or total thrust) can be calculated. This 
use of the pressure probe and the results of the measurements 
are given in Chapter 3. The impact pressure, /O u , can be 
combined with the available velocity profile data to calculate 
the mass flux.you, and kinetic energy flux, . In Chap- 

ter 5, the results of these calculations are presented, a mass, 
momentum and energy balance around the control volume dis- 
cussed and an MPD thruster development program suggested 
based on this new information. 

The investigation described in the following chapters 
lived up to expectations. The new information is significant: 
not only providing an improved understanding of MPD accelerators. 
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but also offering a new method of measuring performance 
which should improve communication between the performance 
and analytical type of investigations. 
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CHAPTER 2 

EXPERIMENTAL APPARATUS 
2 C 1 MPD Accelerator 

The facility used to simulate a pulsed, self field, MPD 
electric thruster is described in detail in Refs. 2-1 and 
2-2. This device accelerates a slug of gas through an MPD 
arc to form a limited duration, but steady flow (after the 
initial transient) exhaust plume which is the subject of the 
present study. A schematic diagram of the facility is shown 
in Fig. 2-1. The elements of this facility which are of 
particular interest are 1) the cylindrical discharge chamber 
(Fig. 2-2), consisting of a 3/4 inch conical tungsten cathode, 
an aluminum anode with a 4 inch orifice and a 2x5 inch ID Pyrex 
chamber wall; 2) the mass injection system (Fig. 2-3), con- 
sisting of the driver and driven sections of a vacuum shock 
tube and three sets of six calibrated orifices for injecting 
argon into the chamber at nominal rates of 1.9, 3.8, 5.9, 23, 
or 36 grams per second; 3) the 160 kilojoule capacitor bank and 
network (Fig. 2-4, Refs. 2-2) for providing a nearly square 
current pulse to the arc; and 4) a 3 ft diameter by 6 ft long 
Plexiglas vacuum tank (Fig. 2-5) for providing a 10 -5 torr 
region into which the exhaust plume may expand without being 
perturbed by EM fields generated by currents in the tank wall. 

The operation sequence of the facility begins with the 
pressurization of the shock tube driver section to 45 or 30 psig 
(depending on the mass flow desired. Ref. 2-1, page 122), 
and charging of the capacitor bank to a voltage from 2 to 6.5 kV, 
depending on the current desired. The shock tube diaphragm is 
ruptured and a portion of the resulting pressure wave is first 
bled into the gas triggered switch and the remainder establishes 
a pressure reservoir at the end of the driven section (Ref. 2-1, 
p. 106) , which in turn feeds the argon gas into the chamber 
through the six calibrated orifices. The length of the bleed 
line to the gas triggered switch is adjusted so that the gas in 
the switch reaches the Paschen voltage, fires, and transfers the 
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capacitor bank potential to the accelerator electrodes only 
after a steady mass flow has been established in the arc 
chamber (Ref. 2-1, pp. 118). A coaxially symmetric discharge 
(Refs. 2-3, 2-4) forms in the discharge chamber which, after an 
initial transient phase, steadily accelerates the incoming gas 
to form an exhaust plume (Fig. 2-6, Ref. 2-3). 

The terminal voltage of the discharge was measured using 
a Tektronix 6013A high voltage probe in conjunction with a 
Tektronix 555 oscilloscope. The probe was calibrated with 
the square wave calibration generator built into the scope. 

Overall accuracy of t 5% is expected (Refs. 2-5, 2-1). 

The discharge current was measured for each shot using a 
Rogowsky coil mounted in the anode barrel (Fig. 2-7) and a 
passive integrator. This coil and integrator were calibrated against 
a precision current transformer using a 555 oscilloscope to record the 
voltages from both the Rogowsky coil and the current transformer. 

The calibration constant was determined to be 1.53 kA/mV, and 
the accuracy estimated to be i 1%. 

2.2 Pressure Transducer Design 

The requirements for a pressure transducer to measure cham- 
ber and jet dynamic pressure in an MPD accelerator are sufficient- 
ly stringent that a special design is necessary. The sensitivity, 
time response, and noise discrimination follow directly from the 
anticipated signal strengths and duration. Pressures on the order 
of ten torr are expected, so a sensitivity of about one torr is 
required, and the noise-equivalent-signal induced in the trans- 
ducer must not exceed one torr to maintain a reasonable signal- 
to-noise ratio. Since the pressure is expected to attain a steady 
value in a few tens of ^isec, and retain it for about one msec, 
it is possible to use an acoustically simple system, where re- 
flected stress waves are held away from the piezoelectric crystal 
in an acoustic delay line until after the time of interest. 

The noise level in a pulsed MPD accelerator is very high, and the 
noise rejection ratio required was found experimentally to be far 
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beyond the capability of commercially available transducer 
systems, so a pressure transducer had to be designed specifi- 
cally to meet these requirements. 

The sensitivity requirement vas easiest to satisfy by 
utilizing a 3/8 inch diameter, well aged, PZT-5, piezoelectric 
crystal. Since PZT-5 has a piezoelectric constant (d^) of 
about 374 picocoulombs/newton (Ref. 2-6), the crystal output 
for a given pressure is about: 

Q = 374x10 s 2.68 x 10 ^ coulombs/N/m^ 

2*200,000 electrons/N/m 2 

A Kistler Model 565 charge amplifier was used to integrate the 
piezoelectric output. Using a 500 picofarad feedback capacitor, 
the output (V) of the charge amplifier fed by the above crystal 
is : 

V 5.36 x 10" 5 V/N/m 2 

With this sensitivity, a pressure of one torr (133 N/m 2 ) re- 
sults in a signal of about 7.13 millivolts, which can be de- 
tected and recorded using a 1A7A preamplifier in a Tektronix 
555 oscilloscope. This crystal, charge amplifier and oscillo- 
scope easily satisfy the requirement that one torr be detect- 
able . 

The noise problem was considerably more difficult to re- 
solve. A study of the noise induced in pressure transducers 
in the pulsed MPD arc environment was carried out using a Kistler 
Model 601 and a very fast transducer (Ref. 2-7) available in the 
lab. Even with these transducers electrically isolated from the 
plasma, the noise-equivalent signal on a one msec time scale 
exceeded one atmosphere, compared with the one torr desired. 

Four major sources of noise were identified. These were: 

1. induced EMF from the rapidly changing magnetic fields? 

2. capacitively induced charge on one or both of the two 
electrical leads from the crystal to the charge amplifier? 
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3. heating of the piezoelectric crystal or its support 
elements, causing thermally induced strains; 

4. and spurious stress waves, principally resulting 
from internal reflections but also from inadequate 
acoustical isolation in the transducer mounting. 

The induced EMF expected was of the order of 10 volts 
per centimeter squared. Since the capacitance of the trans- 
ducer and leads was about 100 picofarads, the induced charge 
from this source will be on the order of 

Q a# 10 V/cm^ X 100 pc/V 04 10^ pc/cm^ 

Since the sensitivity to pressure is about 3.6 pC/torr, to keep 
this source of noise below the equivalent of one torr, the net 
area linked by the changing magnetic field must be kept to less 
than 0.3 mm . With reasonable care in keeping the electrical 
leads close together and using coaxial geometry, this can readily 
be achieved, and the noise from this source in the final design 
was indetfcctable . 

The magnitude of the electrostatic charge induced across 
the crystal was not calculated because of the complex geometry. 
Experimentally, however, it was found very easy to induce a 
charge of more than 3.6 picocoulombs . After considerable trial 
and error, it was found that by forming the outer (ground) con- 
ductor into a completely closed electrostatic shield around the 
crystal and the inner conductor, extending this shield back to 
the charge amplifier, and separating the shield by at least 1/8 
inch from the plasma by potting epoxy, the capacitively induced 
noise from the MPD arc was reduced to a tolerable level. During 
use of the probe, abrupt appearance of noise usually was caused 
by failure of one of these two elements. Either the potting 
epoxy failed and allowed the plasma to get too near the shield, 
or the shield developed a hole. 

Thermally induced strain due to heating of the piezoelectric 
crystal and its support elements proved to be one of the most 
difficult noise sources to remove. Although heat transfer is 
normally characterized by longer time constants, the combination 
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of high crystal sensitivity to thermally induced strain (about 
1000 pC/°K) and the very high plasma temperatures, caused heating 
to be the major source of noise in the Kistler 601 and York (Ref. 
2-7) transducers. The problem was eventually reduced by using a 
variety of techniques reported in the literature, (Ref. 2-10) . 

The principal steps were to change the design from the Kistler 
or York models to one involving two rods, each bonded to one 
face of the crystal, and the crystal edges protected by a thick 
layer of very elastic material. This design isolated the crystal 
from the hot plasma until after the time of interest, and used 
no support structures whose heating and consequent expansion 
would apply stresses to the crystal. It was found later that 
it was necessary also to coat the surfaces of the rods and 
shield the elastic material against radiative heating. Other- 
wise, the energetic radiation generated by the discharge pene- 
trated the transparent rod and potting materials and rapidly 
heated the crystal. The steps taken to reduce this radiative 
heating source of noise were not always successful. Pinholes 
in the paint coating, reflected radiation and erosion of the 
opaque coating allowed some radiation to penetrate, especially 
at the higher current levels. Consequently in the chamber pres- 
sure measurements and the dynamic pressure measurements near the 
discharge, the noise equivalent signal criteria of one torr was 
not always satisfied. 

The two types of acoustic noise were handled separately. 
Acoustic isolation of the probe from the shock tube, building and 
vacuum pump noise was achieved by suspending the probe assembly 
from rubber bands. The pendulum frequency ( Vg/i ) of 40 jisec 
and the tension frequency (\^c/m' ) of 20 |isec achieved were suf- 
ficiently low to filter out that environmental noise with a fre- 
quency comparable to the signal. The low frequency noise which 
passed through the acoustic filters appeared as a baseline shift 
which was a nuisance, but tolerable. This degree of isolation 
from the environment made positioning of the probe difficult, but 
it was found experimentally to be necessary. "O" ring mountings 
on the transducer rods, which would have been much more convenient, 
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were found to produce a noise-equivalent-signal of hundreds of 
torr . 

Internal reflections were minimized by making the probe 
assembly as near a single, homogenous rod as possible. Com- 
plete acoustical homogeniety was impossible since the rod 
material (Plexiglas), electrical contacts and shielding (copper), 
piezoelectric crystal {PZT-5) and bonding material (epoxy), all 
had different acoustic impedances. Similarly the assembly could 
not be a simple rod since some support for the crystal was 
necessary, and because of the thermal and electrostatic shield- 
ing required on the sides of the crystal. By making the non-rod 
materials as thin as possible, however, and by taking considerable 
care to exclude voids which would produce very drastic acoustical 
impedance mismatches in the bonded surfaces, the assembly was 

made effectively homogenous for excitation stress frequencies 
5 

below 10 Hz. Similarly, by using a crystal of the same diameter 
as the rods, using supporting and thermal shielding materials 
with a very low Youngs modulus, and designing the shield to be 
very weak, most of the stress was confined to the rod and crystal. 
To a considerable extent, internal reflections and complex modes 
were thus avoided and the probe acted as a simple homogenous rod, 
with only the longitudinal mode excited by the pressure signals. 

As mentioned earlier, the pressure signals, after a few tens 
of psec of transient behavior, were expected to maintain a steady 
value for the remainder of the current pulse. To follow this 
expected signal, the probe assembly was designed to have a rise 
time of less than 0.1 msec by mounting the crystal 2.0 inches 
from one end of a 38 inch rod assembly. The stress wave generated 
by the pressure pulse on the end of the 2 inch stub would require 
25 psec to traverse the 2 inches of Plexiglas rod (acoustic 
velocity about 2100 m/sec). Any reflections at either crystal 
surface would have a second chance to enter the crystal after 
traversing the stub twice in 50 /isec, so the rise time, even with 
poor coupling between the rod and crystal should be less than 100 
jisec. The remaining 36 inch length of Plexiglas rod on the other 
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face of the crystal provides an acoustic delay line. After 
the stress wave leaves the crystal, it travels to the free end, 
is reflected and travels back to reenter the crystal only after 
about one msec. During the one msec of interest, therefore, 
the crystal is completely free of reflected stresses and in- 
terpretation of the output signal in terms of pressure on the 
stub end is straightforward. For times longer than one msec, 
the output is not simply related to the pressure pulse, but 
this is acceptable since the maximum duration of the MPD cur- 
rent pulse, and consequently any signal of interest is one msec. 

Figure 2-8 shows a schematic of the probe assembly con- 
structed according to the above design criteria. The piezso- 
electric crystal and shield system is emphasized at the expense 
of the simpler but much larger Plexiglas rods. Each element 
in the assembly is identified by an item number which refers 
to Table 1. Table 1 explains the role of each element in de- 
termining the sensitivity, noise and response time of the probe 
assembly . 


2.3 Probe Calibration 

To calibrate the pressure probes it is necessary not only 
to determine the voltage output for a given pressure, but also 
to make sure that the calibration constant did not depend on the 
frequency or amplitude of the pressure over the anticipated 
range. This was accomplished by using four independant cali- 
bration techniques and studying the acoustic properties with two 
excitation techniques. This extensive testing and the many problems 
identified and solved using these tests during the probe develop- 
ment, served to generate considerable confidence that the cali- 
bration constant was valid for the measurements in the plasma. 

2.3.1 Acoustic Properties 

In the acoustic tests, stress waves were excited in the probe 
by two devices; a pendulum and a shock wave. The pendulum con- 
sisted of an 1/8 inch ball bearing suspended from a 6 inch long 
thread. In the test, the ball was held about one centimeter from 
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TABLE I 

Pressure Transducer Construction 

1. 3/8 x 2" Plexiglas stub - This provides thermal and 
electrical separation of the crystal from the plasma: 

It is sufficiently long to remove the crystal from 
the region of strong magnetic fields. It is acous- 
tically simple - only longitudinal stress waves are 
excited. Its acoustic impedance matches the backing 
rod. 

2. Silicone vacuum grease - This is used to achieve a 
good (no void) acoustical coupling between the rod 
and the foil. It is the weakest link in tension, so 
accidental damage will cause failure here, but the 
(no void) acoustical match is easily reformed. 

3o 0.001" annealed and rolled copper foil - This pro- 
vides one of the conductors for carrying the induced 
charge to the charge amplifiers. It also forms part 
of the gaussian electrostatic shield around the crystal and 
center conductor. 

4. GC Electronics #21-1 silver conducting paint - This 
provides electrical contact between the tinned surface 
electrode on the crystal and the Cu Foil. It also 
provides a no void, acoustical coupling. 

5. Piezoelectric (ceramic) crystal, Clevite PZT-5 6050-5 

When strained, this crystal induces electric charges 
on its faces. 

6. 0.001" annealed and rolled copper foil - This 3/8 dia. 
disc conducts the piezo-current from the tinned face 

of the crystal to the center wire of the coaxial output 
wire 0 


Item numbers refer to Fig. 2-8 
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TABLE I (Cont'd) 

7. One mil mylar scotch brand electrical tape - This 
serves to insulate the #6 copper conductor from the 
shield (#9) . The one sticky surface wets the #6 
disc and provides a no void acoustical coupling. 

8. Duco 50/50 formula Epoxe - This provides a no void 
acoustical coupling between the non-sticky side of 
the mylar tape (#7) and the shield (#9). 

9. 0.001 inch annealed and rolled copper foil - This 
foil, together with #3 forms a Gaussian electrostatic 
shield completely enclosing the crystal and center 
conductor. The shield is formed in such a way that 
stresses tending to compress the crystal will bend 
the thin foil. This is to insure that almost all of 
the stress applied deforms and generates a signal 
from the piezoelectric crystal. The #9 foil is 
crimped and soldered to the #3 foil all the way around 
the periphery (except where the center lead penetrates) 
to form a completely closed shield. It is extended 
into a foil shield closely wrapped around the Formvar 
covered center conductor. This shield is then con- 
nected all the way around to the outside conductor of 

a length of RG/58U coaxial cable. At all points the 
separation between the two conductors is minimized to 
reduce the area involved in dB/dt noise pick-up. 

10. Duco 50/50 formula Epoxe - This provides a no-void 
acoustic coupling between the shield and the acoustic 
delay line (#11) and bonds the piezoelectric wafer to 
the rod. 
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TABLE 1 (Cont'd) 

11. 3/8 dia x 36” Plexiglas rod - This rod accepts the stress 
wave propagating through the crystal. The stress propa- 
gates through the rod in the longitudinal mode, is re- 
flected from the free end, and approximately one msec 
later reenters the crystal. This rod, therefore, acts 

as a delay line, keeping reflections from the crystal 
for the time of interest. By keeping the wafer (#4 thru 
#10) thin, and maintaining good acoustical coupling 
throughout the wafer, the #11 and #1 rods are a simple 
acoustic system for wavelengths larger than the wafer 
thickness, and the acoustic impedance mismatch between 
the rods and the various materials in the wafer is un- 
important. 

12. Emerson and Cumings, Inc. Eccosil 2CN potting epoxy - 
This material is used to separate the charged plasma from 
the shield. The large separation reduces the capacitively 
induced shift in the shield (''ground"') potential. It also 
prevents the conductive plasma from discharging to ground 
through the shield. 

13. AWG #30 Formvar insulated copper wire - This wire provides 
the center conductor and insulation of a coaxial cable 
leading to regions of low B field, where conventional 
RG-58/U cable can be used. This cable is made as thin 

as possible both to present as little area between con- 
ductors as possible to reduce dB/dt noise, and to prevent 
the cable from supporting any of the stresses (which should 
be confined to the probe to maintain acoustic simplicity) . 

14. Charge amplifier, Kistler Model 565, S/N 148 with a 550 

pF feed-back capacitor - This amplifier provides an input 

12 

impedance of about 10 ohms and produces a low impedance 
output voltage proportional to the input charge. The 
electronics time constant is about 7 /isec. 
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TABLE I (Cont'd) 


15. Tektronix Model 555 (S/N 00377) dual beam osciloscope 
with type 21 and 21A time base and type L and type 1A7A 
amplifier plug-ins, and Hewlett Packard Model 196A 
Polaroid scope camera - This system was used to display 
and record the voltage (pressure) signals as a function 
of time. 

16. GC Electronics #21-1 silver conducting paint - This was 
used to prevent the intense light generated in the plasma 
from passing through the transparent stub (#1) and heating 
the piezoelectric crystal (#5) . 
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the end of the probe with the string supported from a point 
directly above the probe face, and then released. The ball 
would strike the center of the probe face, and rebound after 
a time (about 50 jisec) determined by the elasticity of the 
impact and the velocity of the stress waves in the ball. The 
natural pendulum frequency of the ball was several seconds, 
so after rebounding the ball remained away from the probe for 
the duration of the test. The result of this impact was to 
induce an impulsive stress wave in the probe. 

Figure 2-9 shows the voltage output from the piezoelectric 
crystal resulting from this impulsive stress wave. The dotted 
line shows the ideal response from a perfectly square stress 
wave into a perfectly elastic uniform, homogenous rod. The dif- 
ference between the two traces is a very sensitive indication of 
the limitations of the probe. 

In Fig. 2-9, the curved initial peak is caused by the small 
contact area and the consequent inelastic and elastic deforma- 
tions. Using a Plexiglas striker, the peak became much broader, 
while using a brass cap on the end of the Plexiglas rod produced 
a much sharper peak. The observed difference between the ideal 
and real results shows a limitation on the test technique in es- 
tablishing the probe response time and in using the momentum 
transferred for an absolute calibration. The deformation of con- 
tact surfaces, however, is not important in measuring gas pres- 
sure which produces a much lower stress per unit area. 

The second small peak shown in Fig. 2-9 is a result of re- 
flection of the stress wave from the Plexiglas - wafer inter- 
faces (Items 2 to 10 in Fig. 2-8). This peak shows the rod is 
not acoustically homogeneous. The ratio of the amplitudes of these 
two peaks is only about 10:1, however, indicating that 90% of the 
stress wave passes through the crystal into the acoustic delay rod 
on the first try. In terms of probe performance this shows good 
response of the probe to rapid pressure changes. 

The third indication of interest in Fig. 2-9 is the long flat 
region from 200 p sec to one msec. This shows that by this time 
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all detectable stress has propagated into the acoustic delay 
rod, and also that this rod is not yet generating reflections. 
(When an "0" ring support or curved delay rod was used, this 
region of the trace was full of wiggles.) With a one msec 
square wave input, therefore, the probe would be free of 
initial transients during this time, and the output voltage 
a simple indication of the stress level. 

The next major feature at about one msec is the returning 
stress wave echo from the end of the delay rod. Since the end 
was free, the compression wave (positive voltage) is reflected 
as a tension wave (negative deflection). This tension wave, 
as shown in Figure 2-9, has been broadened and attenuated by 
the non-linear behavior of irregularities at the ends and all 
surfaces, so before it reaches a maximum, the leading edge has 
propagated into the stub, been reflected from the free end of 
the stub as a compression wave which enters the crystal and 
cancels out the trailing edge of the tension wave. The sum 
of the areas under these reflected waves is about 80% of the 
area of the first spike, indicating only 20% of the momentum 
has been spread out over the rest of the probe by the two free 
end reflections and three passages through the crystal and 
shield. This dissipation is small and indicates the rod is 
acting essentially as a simple homogeneous rod. 

Unfortunately the probe was damaged about half way through 
the experimental measurements and the performance degraded. 
Figure 2-10 shows the results of an echo test made at the com- 
pletion of the experiments. The increase in internal reflec- 
tions is clearly evident from a comparison of Figs. 2-9 and 
2-10. The time from 200 jisec to one msec is not as free from 
disturbances from the initial transient. A comparison of the 
initial and reflected momentum (areas under the curves) showed, 
however, that dissipation was not appreciably effected by the 
increased internal reflections, so the probe voltage was still 
an excellent measure of instantaneous stress. The principal 
damage was to the sensitivity to rapid changes- in pressure. 
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Since the pendulum test could not identify the rise time 
of the probe, a free shock in air was used to measure this 
parameter. Figure 2-11 shows the output of the probe subjected 
to a fast (few jisec) rise in pressure which then decayed in about 
100 /isec to zero. As shown in Figure 2-11, the time required 
for the stress indication to reach its final value is only about 
7 jusec. Since this delay can be entirely ascribed to the elec- 
tronics in the charge amplifier, the shock test shows that the 
response time of the probe is far better than the basic require- 
ment of about 100 fisec. This means that changes in plasma pres- 

5 

sure level at a rate of 10 Hz will be indicated with the same 
calibration constant as slower rates of pressure changes. Fig. 2-12 
shows the decay in rise time caused by the probe damage. 

2.3.2 Calibration Constant 

Four independent methods were used to determine the voltage 
output for a given pressure. Two of these, one based on the 
published piezoelectric constant and the other on momentum 
transferred by collision with a pendulum were relatively in- 
accurate but served to insure that no gross systematic error 
occurred. A third method, using the calculated pressure rise in 
a shock tube was the least subject to systematic error, but also 
the least precise. The fourth method, using a yanked weight was 
the simplest and most precise and its accuracy, or freedom from 
systematic error, was confirmed by the other methods. 

In the yanked weight method, the probe was loaded axially 
with a known weight, and the voltage output measured as the weight 
was rapidly removed. Figure 2-13 shows the voltage history obtained 
by balancing a 200 gm weight on the probe stub end with the other 
end resting on the floor. The scope preamplifier and charge ampli- 
fier were set on D.C. to record accurately the relatively slow 
changes in voltage. The repeatability (which limited the pre- 
cision) of the measurement was good, as indicated in Fig. 2-13, 
but depended on holding the probe perfectly vertical and accurate- 
ly centering the weight directly over the axis of the probe to in- 
sure all of the force acted to stress the piezoelectric crystal 
in compression and did not bend the rod. The accuracy of the cali- 
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YANKED WEIGHT CALIBRATION RESPONSES 
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bration using this technique depends on the assumption that 
the response of the probe to the relatively slow (several msec) 
changes in stress caused by removing the weight, accurately re- 
flects the response to the more rapid pressure changes to be 
measured in the plasma. This assumption was verified with the 
shock tube calibration. 

The yanked weight method proved to be sufficiently simple 
so that it was used each time the probe was removed from the 
vacuum tank, to recheck the calibration constant. The only 
change discovered was after the accidental damage previously 
mentioned, when the calibration constant increased by 19% from 
1.68 to 2.0 volts/200 g weight. This increase occurred because 
the broken potting epoxy (Item 12 in Fig„ 2-8) would no longer 
support any stress. 

Figure 2-14 summarizes the results for a series of tests 
(prior to taking the data) using 5,10,20,50,100 and 200 g 
weights . The bars indicate the range of voltage changes ob- 
tained with each weight and indicate the precision of the 
measurement. The deviations of the measurements from a straight 
line shown on Fig. 2-14 were ascribed to experimental error 
rather than to a non-linear probe response. Neglecting the 
widely scattered 5 g results, and applying the statistical 
test to a least squares fit of the data, the results summar- 
ized in Fig. 2-14 gave a 95% confidence that the probe response 

is within the range of 8.5 t 0 o 2 millivolts per gram weight. 

2 

Since the area of the probe is 0.726 cm , the calibration con- 
stant is 

K = (8.5 t 0.2 X 10 -3 ) (0.726 x 10" 4 )/ 9.8 X 10 -3 
= 6.30 + 0.15 X 10" 5 volts / N/m 2 

= 6.30 i 0.15 volts/ atm 

After the accidental probe damage, a similar series of calibra- 
tions gave K = 7.50 i 0 o 20 volts/ atm. 

In the shock tube calibration method, three pressure probes 
were mounted in the driven section of a 1.2 inch shock tube (Ref. 
2-8) . One probe in the end wall and another 20 cm upstream were 
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used to measure the shock wave velocity. The pressure probe to 
be calibrated was mounted flush with the tube wall, 10 cm from 
the endwall. The shock tube was operated with driven section 
pressures of 1 , ^ 5 ,^, 1 / 10 , 1/20 and 1/100 atm of argon and the 
driver section at zero to 30 psig. 

Figure 2-15 shows a typical data photograph. The bottom 
trace records the voltage output of the probe to be calibrated 
on a 100 /isec/div time scale, and the top trace records the vol- 
tage sum of the outputs of the timing probes on a 20 /asec/div 
time scale delayed exactly 200 psec from excitation of the up- 
stream probe. The apparent ending of the upper trace is caused 
by the shock arriving at the end wall and driving the summed 
voltage off scale. The upper base line contains 10 /isec per 
pulse timing markers to allow accurate shock time of flight de- 
termination. 

Figure 2-16 summarizes the results of the shock tube cali- 
bration. The observed voltages from the data photograph are 
translated into pressure change using the yanked weight cali- 
bration of 6.3 v/atm and plotted against the Mach number squared 
of the shock calculated from the measured time for the shock to 
travel the known distance (20 cm), and using the room temperature 
(15°C) velocity of sound in argon. The bars represent the range 
of voltages observed on a single photograph (see Fig. 2-15) and 
the estimated possible error in reading the driven section pres- 
sure on a mercury manometer. Also plotted on Fig. 2-16 is the 
pressure ratio predicted by the perfect gas relationship 

(m 2 - 1) (Ref. 2-9) 

The triple line indicates the uncertainty in the proper tempera- 
ture (1 5°C) , distance over which the time of flight is measured 
(+ 2 mm) and the molecular weight of air (molecular weight uncer- 
tainty of .05) in the driven section. 

The excellent agreement between the predictions and measure- 
ments of the pressure rise behind a shock based on the yanked 
weight calibration constant proves that the probe response to 
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slow (several msec) and fast (10 ^isec) pressure changes are 
nearly identical. Figure 2-15 shows the easily obtained, high 
precision yanked weight calibration constant is also accurate 
at least within the precision of the shock tube measurements. 

The final two calibration techniques were relatively sub- 
ject to undetermined systematic error, but were useful to re- 
duce the possibility of a major systematic error common to the 
yanked weight and shock tube calibrations. The first of these 

is based on the published average value of the piezoelectric 

— 12 

constant of PZT-5 of 374 X 10 coulombs/newton (Ref. 2-6) . 

The Kistler charge amplifier was found to produce 2 mV/pC, so 
the expected average calibration constant is about 5.4 V/atm 
compared with the 6.3 v/atm for this particular crystal. This 
(15%) is well within reported (Ref. 2-6) variations in this 
material. 

Finally, both the probe and another Plexiglas rod with equal 
mass, m = 80 g, and the same acoustic lengths were coaxially 
suspended by 2.35 meter long strings (1) with their ends just 
touching. The rod was drawn back 5,10,15 and 20 mm (d) and re- 
leased. Figure 2-17 shows the voltage output history from the 
probe. The shape of the rise and decay times indicate that the 
collision was not perfectly elastic and the rod was observed to 
rebound only about one-fifth of the initial separation distance 
(d) . Using these data, the estimated momentum induced in the 
probe (1.2 X ml V 2g/ d ) was compared with the momentum in the 
probe calculated from the area under the voltage history curves. 
(This area in volt-seconds times the ratio of the crystal area, 
0.726 X 10 m , to the yanked weight derived calibration con- 
stant, 6.3 X 10 ^ V/N/m^, gives the momentum of the probe.) 

The ratios of these two momenta were calculated for the 5,10,15 
and 20 mm initial separations as, respectively, 1.02,0.97, 1.00, 
and 0.96 all of which are much closer to one than expected from 
the rough estimate of the rebound distance. This calibration 
also supports the contention that the yanked calibration con- 
stant is indeed a constant. 
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The conclusion is thus that the probe calibration is 
6.3 V/atm, (7.50 after the accidental damage) , independent 
of pressure amplitude or frequency over the range of interest, 
and these values were used with confidence to reduce the plasma 
pressure measurements. 
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CHAPTER 3 

EXHAUST PLUME MOMENTUM FLUX 

The momentum flux in the exhaust plume of the MPD acceler- 
ator was measured with an impact pressure probe. The interaction 
of this probe with the plasma was similar to a pitot tube, so 
the standard pitot corrections were applied to reduce the impact 
pressure data to axial momentum flux. Radial profiles of the jet 
were obtained for eight different combinations of mass flow rate 
and arc current. In addition, centerline momentum flux was 
measured at 1.9, 3.8, 5.9, 23, and 36 grams per second nominal 
mass flow rates for arc currents from 8 to 50 kiloampBres. At 
the "standard" (m = 5.9 g/sec ,1 = 15.6 KA) condition, radial pro- 
files were obtained at 11, 5, and 1 inch from the anode face. All 
profiles were integrated assuming axial symmetry to obtain the 
total momentum content of the exhaust jet. These measurements 
were used to determine the effect of mass flow rate and arc cur- 
rent on the shape and time evolution of the jet and on the thrust 
of the MPD accelerator. 


3.1 Apparatus 

Figure 3-1 shows the apparatus used to measure momentum flux. 
The impact pressure probe behaves essentially as a one inch diam- 
eter, shrouded pitot tube using the pressure probe described in 
the previous chapter to sense the impact pressure. Unusually fast 
response times were achieved by mounting the pressure probe sensing 
surface 1/8 inch behind the probe orifice. 

The 2-3/4 inch Plexiglas tube shown in Fig. 3-rl is used to 
protect and support the pressure sensor and pressure shield. The 
tube prevents buffeting of the sensor by the turbulent jet and de- 
lays any other spurious pressure signals, for example, the rise in 
static pressure on the end of the acoustic delay rod. Consequently 
together with the rubber band pendulum support, the tube acoustical- 
ly isolated the pressure sensor from all acoustic excitation except 
the desired signal for the one msec test. In addition, aluminum 
foil and silver paint on the surface of the tube help shield the 
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piezoelectric crystal from radiant heating by the arc. The 
tube also supports the pressure shield and the rod for ro- 
tating this shield in or out of position from outside the 
sealed vacuum tank. This shield, (see Fig. 3-1) consists 
of a flat, one inch square plate. The place can be po- 
sitioned well away from the pitot tube where it will not 
interfere with the plasma in front of the probe orifice, or 
alternatively, it can be rotated to tightly seal the probe 
orifice. In the latter position, by preventing the plasma 
from reaching the pressure sensor, it allows the noise ampli- 
tude at a given position to be determined immediately after 
a data shot by rotating the shield into position, and then 
taking a duplicate shot. 

The 2-3/4 tube is mounted on a movable platform which is 
suspended from the tank walls. With the tank evacuated, the 
platform axial and radial position may be controlled from out- 
side the tank. Since the tube is nearly as large as the tank, 
the achievable positions of the pitot orifice are limited to 
20 inches from the anode and 10 inches from the centerline. 
Also an obstruction on the platform prevents probing to the 
right of the centerline. 

The pitot pressure history is recorded by a Tektronix 555 
oscilloscope fitted with a Hewlett Packard Model 196A camera 
using Type 47 Polaroid film. The 1A7A pre-amplifier in the 
scope is fed by a Kistler Model 365 charge amplifier with a 
500 pF feedback capacitor. The charge amplifier received the 
pressure signal through an 8 foot length of RG 58/U coaxial 
cable which passes through a vacuum fitting on the tank, 
through Tygon tubing, to connect to a vacuum seal on the pres- 
sure probe coaxial lead. Arc current is detected with the 
Rogowski coil described in Chapter 2 and recorded (on the same 
photograph as the impact pressure history) by the upper beam 
of the 555 oscilloscope. 

3.2 Data Analysis 

The raw data were obtained in the form illustrated in Fig. 
3-2 which shows the arc current, impact pressure and noise 
histories. Each oscilloscope photograph is a record of several 
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shots taken at the same radial and axial position and with the 
same operating conditions of arc current and mass flow rate. 

From this data the axial component of the momentum flux during 
the quasi-steady phase was calculated and the accuracy and pre- 
cision estimated from the multiple shot records for each of the 
positions and set of operating conditions investigated. Then, 
for each axial position, the total momentum flow (thrust) was 
computed by integrating over the jet area, assuming axial 
symmetry. 

The calculation of the momentum flux from the raw data 
involved consideration of the calibration factors, noise sub- 
traction, time selection, conversion of the measured stagnation 
pressures to momentum flux by a pitot tube model, and statis- 
tical analysis of the data scatter. Similarly, the total thrust 
integration involved corrections from the apparent probe position 
to the actual position of the measured impact pressure. The de- 
tails and reasoning for the data reduction procedures are given 
in Appendix 3A. Table II illustrate® these procedures by reducing 
the raw data of Fig. 3-2 (m = 5.9 g/sec, Z = 11 inches and 
R = 0 and 8 inches) to momentum flux at the corrected positions 
and form required for total thrust integration. 

3.3 Results 

The impact pressure histories obtained as a function of 
position over a range of currents and mass flow rates, add appre- 
ciably to the information available concerning the time evolution 
of the axisymmetric flow of accelerated plasma. In„ addition, 
the magnitudes during the quasi-steady phase show the dependence 
of total thrust on current and mass flow rate. This, in turn, 
throws some light on which of the many theories best describe the 
behavior of an MPD accelerator of this geometry and mode of opera- 
tion. 


3.3.1 Operating Conditions 

The matrix of Table III shows the mass flow rates (m) and 
arc currents (I) for which full momentum flux profiles were obtained 
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TABLE II 

Data Reduction Example 

Correction Example* 

1. Time selected: t = 100 jusec for I, 500 ^isec for P 

2. Noise: P (0) — 165 mV - 20 mV = 145 mV 

o 

P (8) = 10 mV - 6 mV = 4 mV 

o 

3. Scaling: I = (1.53 kA/mV) (10.5 mV) = 16 kA 

P Q (0) = (0.159 atm/V) (0.145V) = 0.023 atm 

P q ( 8) = (0.159 atm/V) (0.004V) = 0.00064 atm 

4. Pitot Shock Wave: P(0) = 0.023/0.9 = 0.026 atm 

P (8) = 0.00064/0.9 = 0.0007 atm 

5. Pitot Yaw: P(8) =(0.0007) (0.98) = 0.0007 atm 

6. Pitot Re: P(8) = 0.0006 t 0.0001 atm 

7. Axial Component: P (8)= (0.0006) (cos 8 ) = 0.00048 atm 

z 

8. Axial Flow: P ( 8) = ( 0 . 00048) (cos£) = 0.00039 atm 

zz 

9. Radial Position: P(0) — ► P(0.5") 

P (8) ► P (8. 5" ) 

10. Grad P: P ( 0 . 5) — *»P ( 0 . 5) 

P(8.5)— P(8.36) 


11. 

Symmetry Integral: 

2 JTRP 

zz 

(0.5) = 5 N/inch 



2 itrp„ 
zz 

(8.36) = 1.5 N/inch 

12. 

Precision: 

2 ITrp 

zz 

(0.5) = 5 ± 9 % 



2 fTRP 

zz 

(8.36) = 1.5 ± 115% 



r* 


13. 

Thrust: 

2 ITRP 

dR = 77 N ± 20 % 
zz 


o 


* Raw data from Figure 3-2, m = 5.9 g/sec, Z = 11 inches 



47 


TABLE III 

Operating Conditions 


Profiles at Z = 11 inches 


(g/sec) s \ 

8 

12.3 

15.6 

27 

36 

1.9 

X 



X 


00 

• 

m 


X 




5.9 

X 


X 

X 


23 




X 


36 





X 


Profiles at A = 5.9 g/sec, I = 15.6 kA 

1. Z = 1" 

2. Z = 5" 


Centerline Momentum Flux at Z = 11 inches 

1. For m = 1.9 g/sec, I = 24.7 to 9.4 kA 

2. For m = 3.8 g/sec, I = 30.2 to 7 kA 

3. For m = 5.9 g/sec, I = 30.2 to 11.3 kA 

4. For m = 23 g/sec, I = 41 to 7.1 kA 

5. For m = 36 g/sec, I = 50.5 to 16.5 kA 

The values of the parameters which appear on the diagonal of Table 
III were selected to satisfy the "matched" flow conditions 
(I /m«40 kA .sec/g, see Ref. 3-7). Reference profiles were re- 
peatedly obtained throughout the experimental period for the 
"standard" conditions, I = 15.6 kA, m = 5.9 g/sec. 

In addition to the parameter ranges shown in the matrix, pro- 
files were measured for m = 5.9 g/sec, I = 15.6 kA at 1 and 5 
inches from the anode face to ascertain the jet shape. Finally 
the centerline momentum flux was determined over the operating 
range of arc currents for each of the five available mass flow 
rates. The operating range limits on the arc current were set by 
the limitations of the capacitor power supply system. 


i 
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3.3.2 Jet Description 

Time History A planned delay between the initiation of mass 
and current flows causes the chamber to be filled with an ex- 
panding cloud of neutral gas when the arc is ignited. This 
cloud of gas is ionized and accelerated into the vacuum cham- 
ber (Ref. 3-7) with a relatively slow velocity of less than 
10 4 m/sec (Ref. 3-1) . The leading edge of the accelerated gas, 
therefore, consists of a high density, low velocity region 
which should appear to the impact pressure probe as an initial 
transient. 

Figures 3-2 and 3-3 show the momentum content and radial 
distribution of this initial transient at a distance of 11 
inches from the anode face. In Fig. 3-2a, the initial trans- 
ient is large, indicating most of the swept gas (momentum) is 
accelerated along the axis, but the small initial transient in 
Fig. 3-2b shows some reaches this radius. Since this plasma is 
slower (Ref. 3-8) , but more massive (Ref. 3-7) than the quasi- 
steady jet, it provides an ambient background source from which 
mass may be ingested into the quasi-steady jet. In this way, 
the initial transient provides evidence of a condition which 
might cause the large (100% or more) shot-to-shot variations 
in the momentum flux measured at the edges of the quasi-steady 
jet discussed below. 

After the initial transient, the voltage, current and jet 
velocity (Ref. 3-8) quickly settle to steady values that remain 
constant for at least 0.6 msec. This is the quasi-steady mode 
of operation. Figure 3-2a shows that the momentum flux on the 
centerline is also quasi-steady during this period. (The early 
fluctuations observed in Fig. 3-2a are acoustic oscillations 
in the probe as are the longer lived oscillations in the damaged 
probe as shown in Fig. 3-4 (see Chapter 2)). At the edge of the 
jet, however, the fluctuations are up to 100% of the observed 
pressure. Since the jet Reynolds number as it leaves the anode 
orifice (Sec. 5-4) is about 1.5x10 , these edge fluctuations 
may represent the onset of turbulent interaction of the fast 
jet with the slow ambient plasma. 



VOLTAGE VOLTAGE 


49 



a) CENTERLINE POSITION 



b) 8 INCHES FROM CENTERLINE 


TYPICAL CURRENT , PITOT PRESSURE AND NOISE HISTORIES 


FIGURE 3-2 
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I 

P 

0.0029 atm/DIV 


CURRENT, PITOT AND NOISE HISTORIES 
AT R=6 AND Z =11 in 

FIGURE 3*3 
I 


P 

0.00725 atm/DIV 

SHORT CURRENT PULSE HISTORIES 
FIGURE 3-4 



TIME ( lOOusec/DIV ) 



TIME (lOOasec/DIV) 


AP 25-P-39I 
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The trailing edge of the quasi-steady flow is defined by 
termination of the arc current. Figure 3-4, taken with a 1/2 
msec current pulse, shows that the momentum flux drops to zero 
after the arc is extinguished. The lower velocity associated 
(Ref. 3-8) with the reduced current level over the last 50 
jisec causes this trailing edge to be spread out over 100 psec 
by the time it has traveled the 11 inches to the impact pres- 
sure probe. 

Jet Profile The shape of the quasi-steady jet at 11, 5 and 1 
inches from the chamber is sketched in Fig. 3-5. The main feature 
of interest in this presentation is that the jet leaves the cham- 
ber as a narrow beam about 2 inches in diameter (from a 4 inch 
diameter orifice!) and expands radially into a nearly gaussian 
profile with a half width of nearly six inches at a position 11 
inches from the chamber. Probe interpretation is difficult in 
the very large pressure gradients at Z=5 and Z=1 inches (Section 
3. A. 4), but the best interpretations indicated that the jet edges 
were initially very sharp. When these three profiles were re- 
plotted for the purpose of integrating the thrust, it was found 
that the function 27 Trp (R)/ZR (approximately the derivative 
of the profile) scaled nearly linearly with the Z position. 

The integrated total thrust at all three positions agreed 
within 2%. 

The most rewarding technique for investigating the profile 
shape at Z = 11 inches is shown in Fig. 3-6. To obtain this plot, 
the axial momentum flux measurements were treated as the ordinate 
of a histogram vs. radial position. The fraction of the area of 
the histogram outside a given R was then plotted against R on 
normal probability paper. As is evident from Fig. 3-6, the data 
points fall nearly on a straight line, indicating the distribution 
is nearly Gaussian. The deviations from Gaussian are principally 
at large radii, where the turbulent boundary causes the maximum 
uncertainty in the data. By similarly plotting the profile data 
for all the other mass flow rates and currents, it was found, 
within the precision of the profile measurements, the shape of the 
jet was independent of m and I. The apparent width (standard 
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RADIAL POSITION (INCHES) 


JET PROFILES AT II , 5 a I INCHES 
FROM THE ANODE FACE 


FIGURE 3-5 

AP 25 ■ 4720 
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deviation mO f = 2.3 inches) also remained constant. This provides 

a first order description of the invariant shape of the jet which 

is especially useful in estimating the total thrust of the jet from 

a single impact pressure measurement (P ) on the centerline using 

z z 

the relation: 

Thrust = 2*TO /2 P (R=0) 

zz 

3.3.3 Total Thrust 

Figure 3-7 summarizes the total thrust information obtained 
by integrating the momentum flux profiles. Also shown are the 
nominal mass flow rates for each profile and the total electro- 
magnetic (EM) thrust predicted assuming an axially symmetric cur- 
rent density distribution (Ref. 3-9) . The error bars represent the 
best estimate of the range within which the thrust can be located 
with 95% confidence. These limits are a best estimate rather than 
a precise statistical statement, principally because it was necessary 
to estimate the momentum injected at the edges of the quasi-steady 
jet from the snowplow initial transient. Despite these limits on 
the precision, three significant conclusions can be drawn from the 
summary in Fig. 3-7. 

First and most important, the data show that the total thrust 

scales as the square of the current over the range of parameters 

investigated. This relation has been previously confirmed by thrust 

stand measurements at lower currents and for different geometries 

(Ref. 3-10 and 3-11) . The data reported in Fig. 3-7 confirms that 
2 

T«CI continues to hold up to at least 36 kA for a pulsed MPD 
accelerator with a very large anode orifice (4 inches),; Since 
the theoretical justification of the I relation is based on the 
assumption of axial symmetry, the data also indicate this symmetry 
is maintained ("spokes" are not important). 

Second, the measured total thrust significantly exceeds the 
predicted electromagnetic thrust for all except two of the conditions. 
This result is particularly interesting since the electromagnetic 
thrust component in a real device would be expected to approach the 
theoretical prediction as an upper limit, not exceed the prediction. 

& 
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For example, even a small (few percent) deviation from exact 
axial symmetry would reduce the EM thrust component significant- 
ly. Of the possible real device effects, only a drastic re- 
duction in the cathode attachment radius could increase the EM 
thrust and the experimental evidence (see Ref. 3-2) indicates 
this does not happen in the present accelerator. This leads to 
the conclusion that in this geometry, significant thrust is pro- 
vided by electrothermal acceleration over the entire range of 
parameters. It was expected that at very low currents the MPD 
accelerator would act as an electrothermal arc jet, deriving 
most of its thrust from the expansion of the heated plasma, but 
the indication in Fig. 3-7 that even at the highest currents, 
from 10% to 55% of the total thrust is electrothermal is sur- 
prising. Since this device does not have a material nozzle, 

•A 

this, in turn, indicates the jXB body forces act (favorably for 
a change) to simulate an effective converging - diverging 
nozzle. 


3.3.4 Centerline Momentum Flux 


The momentum flux on the centerline is of particular in- 
terest for two reasons. First it can be measured with much greater 
accuracy (10%) and resolution (5%) than the total thrust, and 
second, it is much easier to determine, requiring only one 
measurement at each condition. The first of these reasons is 
important in that it allows more subtle relationships to be de- 
tected. The second is important because of the apparently in- 
variant shape of the jet (Section 3.3.2). A useful estimate of 
the total thrust can be simply calculated from the single measure- 
ment of the centerline momentum flux. 

Figure 3-8 shows a typical data record for a survey of mo- 
mentum flux on the centerline as a function of arc current, for 
m = 1.9 g/sec and Z = 11 inches from the anode. The first half 
of the traces should be ignored as probe-caused oscillations. 

The latter half of the figure clearly shows the monotonic change 
in momentum flux with arc current (upper set of traces) . 
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Figure 3-9 summarizes the measurements of the centerline 
momentum flux. In this plot, although the general trend is 
identical for all mass flow rates, there is some evidence that 
lower mass flows produce lower centerline momentum flux. The 
effect, if real, is barely within the precision of the present 
measurements . 

Figure 3-9 also shows the centerline momentum to increase 
as the 2.35 power of the arc current rather than the 2.0 power. 
This probably indicates a narrowing of the jet shape at higher 
currents, an effect beyond the precision of the jet shape 
analysis in Section 3.3.2. 



CENTERLINE MOMENTUM FLUX (atm) 
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Appendix 3. A Data Reduction Procedures 

The task of analyzing the data is to convert the oscillo- 
scope photographic records to momentum data as a function of 
position in the jet, mass flow rate and arc current. This task 
divides into four parts. First, the recorded voltage amplitudes 
of the traces are scaled to arc current and impact pressure his- 
tories using the calibration factors reported in Chapter 2. 

Second, the noise is subtracted and a current and pressure 
selected which are representative of the quasi-steady portion of 
the histories, and the precision of the measurements evaluated 
from multiple shot records. Third, the impact pressures are re- 
duced to axial momentum flux values using a standard model of 
supersonic pitot tube behavior. Finally, for a given axial po- 
sition, arc current and mass flow rate, the momentum flux as a 
function of radial position is integrated to obtain the total 
axial momentum in the jet. To illustrate these steps, the data 
recorded in Fig. 3-2 are analyzed below. 

3.A.1 Scaling 

The starting point for data reduction is exemplified in 
Figure 3-2 which shows two typical measurements of the arc cur- 
rent, impact pressure and noise histories. These data were ob- 
tained with a nominal mass flow rate of 5.9 g/sec at a distance of 
11 inches from the anode surface. Figure 3-2a is a record of four 
successive shots with the probe orifice on the nominal centerline 
(R=0) and Figure 3-2b is a record of six successive shots at 8 
inches from that centerline (R=8) . 

The top trace in each photograph is the multiple record of 
the arc current measured relative to the straight baseline (second 
\ trace from the top). The current calibration factor is 7.65 kA/div 
so the currents are about 16 kA. The next trace down is a multiple 
(X 2 in 3-2a, X3 in 3-2b) shot record of the impact pressure with 
calibration factors of, respectively, 0.00795 atm/div and 0.00159 
atm/div. The bottom (light) traces in each figure are the noise 
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obtained by blocking the probe orifice. The baseline shifts 
especially noticeable in Figure 3-2b, are caused by 20 Hz 
vacuum pump acoustic noise. 

3. A. 2 Time, Noise and Precision 

Tirae For both Fig. 3-2a and 3-2b the time near the center 
of the phbtograph (500 jisec after the start of the current pulse) 
was chosen as representative of the quasi— steady impact pressure. 

For earlier times, the initial transient caused by the snowplow 
gas (Ref. 3-1) caused acoustical oscillations in the pressure sen- 
sor (noticable especially in Fig. 3— 2a after t=150 ^isec) super- 
imposed on the EM induced oscillations (noticable especially in 

3— 2b) . For times later than t=800 jisec the returning acoustic 
echo attenuates the strain in the crystal (Chapter 2 ) . The arc 
current amplitude was read at t =100 ^asec, after the initial tran- 
sient but before the RC decay of the passive integrator had a 
chance to appreciably reduce the integrated Rogowsky coil signal. 

Noise To obtain the true impact pressures from the photo- 
graphs, it is necessary to subtract off the noise. For example, 
i * 1 f • 3—2 the average noise on both figures was subtracted from 
the average signal at t=500 psec to obtain average impact pressures dur- 
ing the quasi-steady phase of 0.023 and 0.00064 atm respectively. 

The short (60 psec) initial time with no signal was used to estab- 
lish the baseline level for both noise and signal amplitudes. The 
noise at t=500 )isec is principally from radiant heating of the 
piezoelectric crystal as can be recognized from the steady in- 
crease with time until the arc is extinguished. 

Precision A statistical analysis of the shot— to— shot varia- 
tion was carried out for one profile to establish the precision of 
the measurements. As may be inferred from Figs. 3-2, this precision 
was nearly independent of the radial position. The sample standard 
deviations for 9 signal minus 9 noise shots at R = 0, 1^, 3 , 4, 5 , 

6 , 7, 8 and 9 inches were, respectively 2.02, 3.15, 2.25, 1.13, 1.86, 
0.19, 0.52, 1.05 and 1.05 X 10 ^ atm. Applying the statistical t 
test, the 95% confidence limits for the true mean have a range about 
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the sample mean of, respectively, 1.5, 2.4, 1.7, 0.85, 1.4, 

-3 

0.14, 0.4, 0.8 and 0.8 X 10 atm. Since the impact pressures 

_ 3 

for this profile varied from 23 X 10 on the center 

-3 

0.6 X 10 atm at R=8, the precision varies from 5.8% to 114%. 
Because of this large imprecision at the edge of the jet, the 
total momentum flux in the jet can be determined only within 
about 30% at the 95% confidence level. This is a significant 
limitation on the results that can be expected from the present 
technique. 

3. A. 3 Pitot Corrections 

Conversion of the impact pressures to .momentum flux was based 
on the standard supersonic pitot model (Ref. 3-2) . This model 
assumes a curved shock forms in front of the pitot orifice. The 
shock location and its effect on stagnation pressure is determined 
by the probe diameter, free stream Mach number, transverse velocity 
gradient, probe Reynolds number and yaw angle. The model assumes 
the plasma on the stagnation streamline passes through the 
shock with a normal incidence and is brought to rest isentropically 
in the subsonic region aft of the shock. 

Supersonic Using this model, the pitot or isentropic stagna- 
tion pressure aft of the shock, P Q , is related to the static (thermo- 
dynamic) free stream pressure, p , and the free stream Mach number, 

s 

M by the Rayleigh supersonic pitot tube formula, (Ref. 3-2) , 



Since the total momentum flux in the free stream, P, is related 

2 2 

to the thermodynamic pressure p by P = p +k?U = p(l +y M ), 

S S f s 

the pitot pressure and momentum flux are related by 



Fortunately, this expression is very insensitive to the exact value 
of both 2T and M, since neither quantity is known in the plasma 
jet with great accuracy. For Mach numbers from 2 to « 0 , and with 
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V near 5/3 in and after the shock and slightly lower in the 
free stream, the above relation reduces to 

P = P / 0.9 

with an estimated uncertainty of 5%. 

Reynolds Number Based on the extensive empirical experience 
with pitot tubes, three additional corrections to the impact pres- 
sure were considered. These were viscous losses, yaw and resolu- 
tion of the pitot pressure into axial and radial components. The 
effect of viscosity on pitot pressures has been extensively in- 
vestigated as a function of Reynolds number. The results, sum- 
marized in Refs. 3-3 and 3-4 indicate a negligible correction is 
required for probe Reynolds numbers greater than 20. The estimated 
density, velocity and viscosity near the center of the jet are of 

the order of, respectively, 5 X 10 ^ (Chap. 5), 2.5 X 10 ^(Ref.3-8) 

-5 

and 10 (Ref. 3-6) in MKS units, so the probe Reynolds number is 
above 300. Near the edge of the jet, however, the estimated 
Reynolds number drops to about 14, so the indicated pitot pres- 
sure is too high. The data near the edge of the jet are suffi- 
ciently uncertain because of measurement imprecision, making a 
precise correction factor meaningless. Consequently, an estimated 
viscosity correction was based on the empirical correction factor 
1 + 5/ Re , reported in Refs. 3-4 and 3-6. Applied to the present 
data, this works out that 20 ± 20 N/m z should be subtracted from 
the momentum flux in the outer regions of the jet. 

Yaw The axis of the impact pressure probe was maintained 
parallel to the jet axis, while the velocity vector direction de- 
pends on the radial position. Experience with pitot tubes has 
shown that the recorded stagnation pressure falls off as the yaw 
angle between the probe axis and the velocity vector increases. 

The top curve in Fig. 3-10 shows the predicted (Ref. 3-4) yaw re- 
sponse of a shrouded probe as a function of the yaw angle. The 
barred lines show the results of an attempt to experimentally 
verify this function for our probe by changing the probe angle but 
keeping the orifice location fixed on the jet centerline. The 
broad spread and the limited angular range of the data were a con- 
sequence of the limited space and motional control in the vacuum 
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EFFECT OF YAW ANGLE ( 9 ) ON PROBE 
RESPONSE TO AXIAL MOMENTUM 


FIGURE 3-10 

AP25 • 4721 
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system. Still, these data support the contention adopted that 
the yaw response is given by the upper curve. 

Tensor With the yaw correction, the calculated momentum flux 

is the magnitude of the flux tensor P = p I + /° uu. This is the 

quantity of momentum in the u direction which flows per second 

through an area perpendicular to u. The axial component (z) of 

this flux tensor through an area with a direction a is p la + 

^Ou 2 u-z u*a where the " * " symbol indicates a unit vector. Only 

the axial component was calculated and the integration was carried 

out over an area normal to the jet axis (a = z) , so only the flux 

element P zz = p +/Ou 2 cos 2 0 (where u*z = cos0 ) was of interest. 

Since the static pressure (especially near the jet edges where 

cos 2 ^ is significant) is small, the term (1-cos^d )p was neglected 

and the tensor element P calculated from the relation 

z z 

P = | P I cos^d 

This conversion factor is also plotted as the lower curve in Fig. 
3-10. In practice, the yaw and tensor corrections were calculated 
simultaneously by using the difference between the two curves on 
Fig. 3-3, Y(d ) , times the measured impact pressure. Fortunately, 
the uncertainty in this correction is principally associated with 
©>30° (R>6), since the yaw correction of zero for smaller angles 
has been well established (Ref. 3-4). For R> 6, the estimated un- 
certainty is 20% of the correction factor. 

Applying these transformations to the data of Figs. 3-2, the 
mean, momentum flux on the centerline P (0) = 23 X 10 3 atm/0.9 = 

26 x 10 J atm and the mean momentum flux at R=8, P (8) = (0.6X10 ) 

0 -R -3 

(0.83) /0. 9 minus 20 X 10 = 0.26 X 10 atm. Using the square 

root of the sum of the squares to combine the uncertainties, the 

true mean of P (0) is expected to lie within the range 23.5 to 
~ z z 

28.5 X 10“ J atm (5% scale, 5% supersonic, 5.8% precision, total = 
8.6%). Similarly the true mean P zz (8) is expected to lie within 
the range 0.56 to -0.04 X 10 -3 atm (5% scale, 5% supersonic, 114% 
precision, 4% yaw, 28% Re, total = 115%) 
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3. A, 4 Integrated Momentum 

The final step in data reduction, is to obtain the total 
momentum flow rate or thrust of the jet by integrating the flux 
over the jet cross-section. Axial symmetry was assumed so this 
integral reduces to (2^*RP )dR. In this form, it is apparant 
that the accuracy of the radial position is just as important as 
accuracy of the momentum flux in determining the thrust of the 
jet. 

Two corrections to the radial position values are required. 

The first correction is a consequence of the effect of the radial 
gradient in the momentum flux on the pitot measurement. This 
problem has been encountered frequently in boundary layer measure- 
ments with pitot tubes, and was studied in plasmas in the work re- 
ported in Reference 3-6. In these studies it was found empirical- 
ly that the effect of an impact pressure gradient normal to the 
axis of the pitot tube was to shift the effective orifice center 
a distance d = /D7P/P + d 5TP D/4 toward the region of higher 
momentum flux. (D is the probe diameter and P the pitot pres- 
sure.) This correction, when applied to our data, resulted in a 
shift of at most 0.2 inches, but the effect of this shift on in- 
tegrated momentum was significant since it narrowed the measured 
size of the jet. 

Centerline Correction The second correction was required to 
locate the true centerline of the jet with respect to the nominal 
R=0 position of the platform. This was very difficult to determine 
beforehand using physical measurements, because the apparatus is not 
designed to make the platform motion along the vacuum tank axis, the 
alignment of the chamber and tank axis, and the chamber and jet 
axis alignment all accurate within the 1/4 inch required. Further, 
the importance of the accuracy of this measurement was not appre- 
ciated until after all of the data were collected. Consequently it 
was necessary to use properties of the data to make the necessary 
corrections. 



67 


Figure 3-11 shows the raw data prior to this correction. 

The integrand of the thrust integral, Iff'RP , is plotted against 
the platform position for the six radial profiles obtained at the 
"standard" conditions (5.9 g/sec & 16 kA.)Of particular interest 
are the two sets of data connected with heavy lines in Fig. 3-11. 
One, (data symbol "B") is defined by the mean of 9 shots at each 
position (the series also used to determine reproducibility) . Prior 
to collecting these data, the edges of the jet were approximately 
located and this information used to set the platform R=0 position. 
This set is therefore the most accurately centered data available. 
The other curve of particular interest (data symbol "o") is defined 
by the one series of measurements where special arrangements were 
made to extend the profile to both sides of the nominal centerline. 
If the jet centerline and R=0 coincided in this case, the 2&RP ZZ 
curve should be symmetric about R=0. Instead, as shown on Fig. 
3-11, the maximum of the "o" curve for R< 0 is lower and occurs at 
a smaller R than the standard curve ("B") and the portion of the 
profile for R>0. A study of the behavior of 27/'RG, where G is 
a Gaussian curve, showed that these are the properties expected 
with this ordinate (which is the derivative of the Gaussian) with 
an incorrect R=0. The curve maximum will move down and to smaller 
R or up and to larger R depending on the direction of the error 
in R=0. With these properties in mind, the R=0 for the full pro- 
file curve ("o") was shifted 0.5 inches to equalize the magnitude 
of the maxima on both sides of the centerline. This also caused 
the radial location of tbe maxima to coincide with that of the 
"p" curve at about 2.2 inches. This confirmed that the correc- 
tion was successful and necessary. Next, the remainder of profiles 
were similarly adjusted so that the radial location of their maxima 
were also at 2.2 inches. This required, at most, an adjustment of 
3/4 inches. The corrected data for the five profiles are plotted 
in Figure 3-12. All data in this display lie close to the standard 
data, as they should since all are measurements at the same mass 
flow rate and current. The scatter can be completely explained by 
the uncertainties previously noted. 
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RAW DATA THRUST PROFILE DERIVATIVES 

FIGURE 3-11 



• 

o 

o 


£ 

o 

.#/ \_t • 

/V V □ 

!/? '°^0‘v A 

•Jj x v 

r n ^sN». ^ ? 

1 1 

3-2-1 i 

8 

Z 1 1 1 i 1 1 <Sa-s.,. ° 

01 23456789 10 


ADJUSTED t THRUST PROFILE DERIVATIVES 

FIGURE 3-12 


AP 25 -4717 


69 


Finally, this correction was applied to the profiles at 
other mass flow rates and currents. Fortunately, it was found 
by checking back through the lab notebook, that most of these 
were taken after a standard condition run (" 57 ") that had re- 
quired no correction. The probe had not been disturbed between 
this "57" run and the 8 subsequent runs so no correction was 
applied to the R=0 for these profiles. One of the remaining 
immediately followed a standard condition run "0", so the 
centerline correction was known. The final profile, taken at 
36g/secand38 kA, was immediately followed and preceded by 
removal of the probe, so no guide to the correction required 

was available. The location of the 2 J^RP „ maximum indicated 

zz 

a correction of 1/2 inches, so this correction to the R=0 was 
used and the limits of uncertainty increased correspondingly 
in reporting the thrust. 

Precision An error analysis of the integration of the pro- 
files show two compensating effects. The uncertainty of each 
point used to calculate the thrust is increased by the uncertain- 
ty in R. This is easily appreciated by considering the thrust 
integrand, . In analyzing the data, an uncertainty in R 

of ± 1/8 inch was assumed. On the other hand, the profiles form 
a pattern, so a least squares fit of the profile would reduce the 
uncertainty of the integrated thrust. Rather than actually using 
a least squares fit, the same effect was achieved (with a mathe- 
matical complexity more in keeping with, the rest of the precision 
analysis) by drawing a smooth mean curve through the mean of each 
point. Then the sum of all point uncertainties was added to half 
the sum of the squares of the variations of the mean points from 
the curve, and the square root of this sum divided by the number 
of points in the profile to give an uncertainty in the ordinant of 
the mean curve. The two limiting curves were drawn and graph- 
ically integrated to give two values of thrust, which are expected 
to bracket the value of the true thrust in the jet. The mean and 
limiting curves are shown on Fig. 3-12. 
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CHAPTER 4 
CHAMBER PRESSURE 

The pressure in the arc chamber of the MPD accelerator was 
measured using a static pressure tap in the side wall. Sources 
of systematic error in the measurement were identified and com- 
pensated. Data were obtained for five nominal mass flow rates 
over a range of arc currents from 8 to 50 kA. With the aid of 
several assumptions, the aerodynamic or electrothermal component 
of thrust was calculated from the chamber pressure and the ratio 
of electrothermal to electromagnetic thrusts obtained as a func- 
tion of m and I. 


4 . 1 Apparatus 

Figure 4-1 shows the apparatus used to measure chamber pres- 
sure. The probe described in Chapter 2 is used as the sensor in 
a h inch ID pressure tap in the chamber wall. The hole and its 
surroundings were painted with reflecting paint to reduce the 
amount of radiant heating of the piezoelectric crystal by the arc. 
A rubber band, pendulum support for the probe and the hole clear- 
ance (3/8 inch probe in a 1/2 inch hole) isolated the sensor from 
acoustic noise. Fast response was obtained by mounting the pres- 
sure sensitive face of the probe 1/8 inch from the interior cham- 
ber wall. The residual acoustic, electromagnetic, electrostatic 
and thermal noise was measured by blocking the static pressure 
tap with mylar tape backed with aluminum foil. This also blocked 
the thermal noise caused by arc chamber radiation, but this type 
noise could be separately identified on the total signal responses. 

Arc current and the terminal voltage were measured for each 
shot with the Rogowosky coil and Tektronix voltage probe described 
in Chapter 2. The arc current history was traced by the upper 
beam of the oscilloscope and recorded on the same photograph as 
the pressure history. The terminal voltage history was recorded 
with a second scope and camera. 
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FIGURE 4- 1 

AP25- 4743 



CHAMBER PRESSURE PROBE 
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4.2 Data Analysis 

The data analysis followed essentially the same pattern 
as that used to reduce the impact pressure measurement: 1) scal- 
ing, 2) noise elimination and selecting a representative time, 

3) static pressure tap corrections to get chamber pressure and, 

4) calculating the electrothermal thrust component. Steps 1 
and 2 were essentially identical so just the salient points 
are given below. Steps 3 and 4, however, are discussed in de- 
tail since the procedures for reducing static pressure data 
differ significantly from procedures for reducing impact pres- 
sure data. 

4.2.1 Scaling 

Figure 4-2 shows a typical data record of the chamber pres- 
sure (lower trace) and arc current (upper trace) histories. The 
current scale shown was derived from the Rogowsky coil calibration 
of 1.53 kA/mV and the pressure scale from the pressure probe cali- 
bration of 6.3 V/atm. 

4.2.2 Precision, Noise and Time 

The precision of the chamber pressure measurements was not 
formally determined. Successive shots (Figs. 4-2 and 4-3) showed 
that shot-to-shot variations were sufficiently low (compared to 
other uncertainties) to be neglected. 

The noise was measured by blocking the static pressure tap 
with tape so it was necessary to open the vacuum system between 
signal and noise shots. Consequently, noise was recorded on a 
separate photograph from the chamber pressure data. Figure 4-4 
shows such a noise record corresponding to the chamber pressure 
data shown in Figs. 4-2 and 4-3. The bottom trace on Fig. 4-3 
has been manually scratched on the photograph to show the net 
pressure history, obtained by point by point subtraction of the 
(Fig. 4-4) noise from Fig. 4-3. 

The time selected on each oscilloscope photograph as repre- 
sentative of the quasi-steady phase, was determined by the type 
and amount of noise present. For most photographs, the time 
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FIGURE 4-2 
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selected was 350 jisec after the start of the current pulse. The 
initial 25 psec no-signal portion of the trace established the 
zero pressure base line. 

4.2.3 Static Pressure Corrections 

Normally there is little question that a static pressure tap 
measures the local free stream pressure (Ref. 4-1) . Some experi- 
mental check of this point seemed desirable, however, because of 
the presence of electromagnetic acceleration in the chamber as 
well as the inevitable gas leakage in the annular space between 
the 3/8 inch probe and the 1/2 inch hole. Consequently, repeti- 
tive measurements were taken to determine the effect of blocking 
the annular space with plastic sponge, and moving the probe to 
different positions along the hole axis. The pressure indication 
did not change when the sponge completely blocked leakage around 
the probe. Acoustic noise increased considerably, but changes in 
pressure indication of greater than about 10% would still have been 
detectable. 

Measurements of the chamber pressure also showed no change 
with the face of the pressure sensor located at five different 
positions relative to the interior surface of the chamber. The 
positions were with the probe face: protruding 1/8 inch into the 
chamber? flush with the chamber interior surface, and withdrawn 
1/8" , 2/8" , and 3/8"back into the hole. This insensitivity to changes 
in the chamber boundary indicates the absence of rapid flow or 
strong pressure gradients which would have required further static 
pressure corrections. 

4.2.4 Electrothermal Thrust 

The contribution of gasdynamic forces to the total thrust of 
the MPD accelerator is the vector integral of the normal pressure 
over all solid surfaces plus the vector integral of the viscous 
forces over the same surfaces. Measurement of these forces as a 
function of position is beyond the scope of this work. Instead, 
several drastic but reasonable assumptions were made to provide a 
plausible estimate of the force exerted on the plasma by the sur- 
faces of the arc chamber. 
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The first assumption was that the viscous forces were 

negligible. In MKS units, the hot plasma viscosity is less than 

10 “ ^ (Ref. 4-3), the velocity about 10^ (Ref. 4-2) so with a 1 mm 

_2 2 

boundary layer the viscous force would be less than 10 N/cm . 

2 

Since the surface areas are tens of cm , and the total thrust of 
the order of 100 N, the assumption of negligible viscous loss 
seems reasonable. 

The second assumption was that the pressure on the chamber 
walls was uniform. Some support for this assumption was offered 
by the experiments described in Section 4.2.3 which indicated the 
ionizing current near the wall did not generate rapid flow or high 
pressure gradients. The current pattern in the arc chamber (Ref. 
4-2) indicates any pressure difference along the anode walls due 
to jxB forces would be concentrated near the walls and would have 
been detected by these experiments. The details of the plasma 
flow are not known, but it seems reasonable that flow pressure 
gradients in the volume behind 80% of the current connecting the 
chamber walls will be small since the flow area is large and the 
jxB accelerating forces are small. 

The final assumption is that the force exerted on the 200 
cm front face of the anode is negligible. This assumption is 
justified by the experimental measurements of the plasma density 
and temperature in Ref. 4-4. The maximum density on the face of 
the anode is about 10 /m , so even with an ion temperature of 
5 eV, the force would be less than 2 Newtons, quite small compared 
with the 100 N total thrust. 

4.3 Results 

Figure 4-3 shows that, as in the case of the terminal and 
field measurements, the chamber pressure enters a quasi-steady 
mode after an initial starting transient. The duration of the 
transient agrees well with the duration of the terminal voltage 
and magnetic field transients. 

Figure 4-5 shows a typical data photograph from a survey 
of chamber pressure as a function of mass flow rate and current. 
Each of the lower traces represent a single shot at the current 
shown in the corresponding upper trace. All data on a single 
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SURVEY OF P ch AS A FUNCTION OF ARC CURRENT 


FIGURE 4-5 

AP25- P-400 
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photograph are obtained using a single mass flow rate. Figure 
4-6 summarizes all of this chamber pressure survey data. On this 
figure, each of the data traces has been reduced to a single 
chamber pressure value by the methods of Section 4.2 and is plot- 
ted vs. the arc current (I) using five symbols to differentiate 
between the five different mass flow rates (m) . Immediately 
apparent from Figure 4-6 is that chamber pressure is strongly 
related to arc current and much more weakly dependent on mass 
flow rate. 

To express these relationships empirically, a least squares 

fit to the data for each m was used to find that chamber pres— 

1 54 

sure scaled as (I) * for all mass flow rates. Each chamber 

pressure measurement data point was then divided by I 1 * 54 to 
eliminate the arc current dependance, and the resulting numbers 
plotted vs. m as shown in Fig. 4.7. The dots in this figure 
represent a single datum, and the "X" ' s the mean at each mass 
flow rate condition. This plot shows that the m dependence of 
chamber pressure is weak, but unambigious. Again a least squares 
fit to the data was used to express an empirical relation for 
chamber pressure (P ch ) in MKS units of 

P ch = 3.02 t 0.05 X 10" 3 (I) 1 * 54 (m) 0,294 

Where the precision limits indicate the range of the mean value 
* „ ,_1.54 . 0.294 

of ch' 1 m within 95% confidence limits. Note, how- 

ever, that a least squares fit does not prove the exponents are 
exactly the values given. 

The electrothermal thrust component (T efc ) expected with these 
chamber pressures, is (using the assumptions discussed in Section 
4.2.4) the chamber pressure times the throat, or anode orifice area 
of 81 cm . Using the empirical relation above (and rounding off 
the exponents) this calculation gives an electrothermal thrust (in 
MKS units) 

T et = 2.45 X 10" 5 I 1 * 5 m 0,3 

or numerically, the electrothermal thrust varies from 5 to 100 N 
over the range of I and m investigated. 


CHAMBER PRESSURE (atm) 
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REDUCED CHAMBER PRESSURE 


FIGURE 4-7 
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The above data provides an answer to the long standing ques 
tion concerning the electrothermal fraction of the total thrust 
of the MPD accelerator. For each of the total thrust measure- 
ments reported in Fig. 3—7, this ratio was calculated using the 
corresponding electrothermal thrust from the data on Fig. 4-6. 
This ratio is plotted in Fig. 4—8 as a function of the parameter 
I. Over the entire range of parameters investigated, about 20% 
of the thrust is caused by plasma pressure on the chamber walls. 

A plausible reason for this fraction to remain nearly con- 
stant is given by an extrapolation of the data using the em- 
pirical relation developed above. Assuming (as verified in 
Chapter 5) that the total thrust is the sum of the electrother- 
mal and electromagnetic )i Q I 2 (In r a /r c + 3/4) /4 T) components, 
the fraction (F) of electrothermal to total thrust is given by 

2.45X10" 5 I 1,5 m 0,3 

2.4X10 -7 ! 2 + 2.45X10" 5 ! 1 * 5 A 0,3 


= 1 + 


(I 2 Zml 

100 m 


1/4 

0705 


-1 


This relation is also plotted in Fig. 4-9 and suggests that the 

A 

range of i/m covered by the experiments was not sufficient to 
show large changes in the fraction F. This relation also pre- 
dicts the arc jet-to-MPD transition which was the original reason 
for the interest in the fraction F, 

It is interesting to note the presence of the parameter 
I /A in the equation for F. This parameter also appears prom- 
inently in the best available theory of MPD accelerator operation 
(Ref. 4-7) , which is based on a minimum principal. The theory 
has not yet been developed to the point where a direct compari- 
son with the experiments is possible, but the appearance of the 
I /A is encouraging. 
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ELECTROTHERMAL CONTRIBUTION 


FIGURE 4-8 

AP25 • 4740 



RATIO OF ELECTROTHERMAL TO TOTAL THRUST 
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RATIO OF ELECTROTHERMAL TO TOTAL THRUST 


FIGURE 4-9 

AP25- 4741 
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CHAPTER 5 

MASS, MOMENTUM AND ENERGY FLOW 

The measurements reported in Chapters 3 and 4 combine 

with earlier data to provide a major increase in our knowledge 

of mass, momentum and energy flow from the MPD accelerator, 

2 

Since the momentum flux is /Ou , dividing by the local velocity 
(Ref. 5-1) gives the mass flux profile of the jet. The integral 
over this profile compared with the total inlet mass flow rate 
(Ref. 5-2), shows appreciable mass is bypassing the jet. The 
calculated electromagnetic thrust (Ref. 5-3) plus the measured 
electrothermal thrust (Chap. 4) equals the jet axial momentum 
content (within the accuracy of the data of Chapter 3) indicat- 
ing the major sources and sinks of momentum have been identified 
and measured. The energy of streaming motion is^Ou times ^ 
the velocity, so the spatial distribution of this major fraction 
of the total energy flow is identified, providing a basis for 
improved estimates of the remaining energy sinks. 

The engineering implications of this knowledge of mass, 
momentum and energy flow follow immediately. The mass utiliza- 
tion efficiency (60%) , specific impulse (1400 sec) , thrust ef- 
ficiency (21 Yo) and other engineering parameters describing the 
performance of the accelerator as a space thruster are easily 
calculated. More important, the position and time resolved in- 
formation shows the reason for the particular values of the per- 
formance parameters. This, in turn, suggests what must be done 
to improve performance and provides clear directions for an ex- 
perimental development program. 

5.1 Mass Flow 

Figure 5-1 shows the axial mass flux profiles at 11 and 
5 inches from the anode which are produced by an arc current 
of 15.6 kA and a mass flow rate to the accelerator of 5.9 g/sec.. 
The data reported in this figure were calculated by dividing the 
total momentum flux taken from the smoothed curves discussed in 
Chap. 3 by the velocity taken from the profiles of Fig. 5-2 , (see 
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Ref. 5-l) at each point and multiplying by the cosine of the 
velocity vector. This method for calculating mass flux is 
based on the assumptions: 1) that the Mach number is high, 
so the impact pressure ( p s t a tic + ^° u ' was Primarily /©u ? 

2) that the velocity vector pointed from the center of the 
anode, and 3) that the extrapolation of the velocity profile 
shown as a dotted line in Fig. 5-2 is valid. A quantitative 
calculation of the uncertainty introduced by these assumptions 
proved to be infeasible, so the accuracy of the mass flux re- 
ported in Fig. 5-1 is unknown. A "best judgement" is that 
the mass flux is correct to 20% near the centerline and 50% 
near the edge of the jet since the sensitivity to uncertain- 
ties in momentum and velocity is linear (eg. 10% error in 
velocity implies 10% error in mass flux) . 


Figure 5-3 illustrates one of the most informative ways 
of displaying the mass flow data. The mass flux profile was 
graphically integrated (using 2lfR as a weighting factor as 


y 

y 

y 




FIGURE 5-3 



86 


discussed in Chap. 3) from R=0 to R. The radial positions 
for which this integral was 10% through 60% of the nominal 
inlet mass flow rate are plotted in Fig. 5-3 for the two 
axial positions. Points representing the same percentage 
define the mass streamlines, which are labledin Fig. 5-3 
with the fraction of the (nominal) total mass flowing in the 
conical volume bounded by the streamline. 

The data as displayed in Fig. 5-3 indicate the steady 
state mass utilization is poor. Only 60% of the mass carries 
appreciable momentum and only a central core moves at the 
highest velocity (see Fig. 5-2) . This conclusion, if true, 
is so important to a thruster development program, that it is 
worthwhile to consider the supporting evidence in detail. 

The value for the inlet mass flow rate is the most 
questionable of the data indicating poor mass utilization. 

The nominal value used in the calculations was obtained by a 
choked orifice, mass flow rate calculation (Ref. 5-2) . The 
accuracy of this calculation is unknown, since it was based 
on a reasonable but nevertheless assumed reservoir temperature 
of 230°K and a reservoir pressure measurement that was probably 
subject to considerable acoustic and heating noise (see Chap. 

2) . An attempt to check this cold flow m is illustrated in 
Fig. 5-4. These data were obtained by blocking the anode orifice 
to provide a known limited volume receptacle for the inlet mass 
flow and used the chamber pressure measurement arrangement 
(Fig. 4-1) to measure the time rate of increase in the chamber 
pressure caused by the cold gas flow. The mass flow rate (m) 
in such an arrangement is the time rate of change of density (/O) 
multiplied by the chamber volume (V) so the measured dP/dt and 
the desired m are related by 

dP/dt = RTm/V + yORdT/dt. 

Again the temperature is not known so the mass flow rate calcu- 
lation depends on an assumption. Assuming, however, that the 
chamber temperature is not a function of time and is times the 
assumed shock tube reservoir temperature of 230°K, the data of 
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Fig. 5-4 indicate a mass flow rate of 15 g/sec. This estimate, 
then, suggests the mass utilization is even worse than shown 
on Fig. 5-3. Similarly any corrections due to mass ablation 
from the chamber walls, or injestion by the jet from ambient 
plasma (Chap. 3) will make the calculated mass utilization 
even worse. 

The velocity profiles as shown in Fig. 5-2 are based on 
time-of-f light measurements and also may contain some systematic 
error. If the velocity were lowered as suggested by Doppler 
shift measurements (Ref. 5-4) , the calculated mass in the jet 
would be larger and the mass utilization better. The Doppler 
shift measurements, however, measure mass averaged velocities, 
so the values reported in Ref. 5-4 which are much less than the 
peak velocity measurements reported in Ref. 5-1 tend more to 
support the proposition of poor mass utilization rather than 
dispute the peak velocities of 25 km/sec. In addition, a veloc- 
ity profile sufficiently smaller (1/3) to imply the injected 
mass is in the jet would reduce the calculated streaming energy 
(Sect. 5-3) to an unreasonably low value. The velocity profile 
as used is probably sufficiently accurate to support the conclu- 
sion of poor mass utilization. The accuracy of the momentum 
measurements was discussed in Chap. 2, and a 1/3 larger momentum 
flux, especially near the centerline seems improbable. 

On balance then, the evidence indicates an appreciable 
fraction of the inlet mass flow is not included in the high 
velocity jet. As suggested in Fig. 5-3 by the dotted line, the 
remaining mass probably is ejected from the 1 inch annular space 
between the 2 inch diameter jet and the anode lip. rather than 
retained in the chamber. The alternate, retention of the excess 
mass in the chamber until the conclusion of the current pulse, 
should cause an increase in the chamber pressure which is not 
observed in the quasi-steady phase (Fig. 4-3) . Also the measured 
chamber pressure alone (0.02 atm), without any electromagnetic 
acceleration, would force the missing 2 g/sec in choked flow 
through the 60 cm annular space between the jet and the anode 
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lip. Either this flow exists or there is some balancing 
force for which no evidence has been found. 

5 . 2 Momentum Flow 

Figure 5-5 shows the axial momentum flow pattern for 
a nominal input of 5.9 gm/sec and an arc current of 15.6 kA. 

This pattern was derived from the momentum flux profiles re- 
ported in Chap. 3 in the same way as the mass flow streamlines 
were calculated from the mass flux profiles. That is, the 
weighted integral from R=0 to R of the average flux profiles 
at 11, 5 and 1 inches from the anode were expressed as a frac- 
tion of the total thrust, plotted against R, the radial positions 
for 10% through 90% crossplotted to give the points marked by an 
in Fig. 5-5, and the points connected to define the momentum 
streamlines. 



AXIAL MOMENTUM FLOW STREAMLINES 


FIGURE 5-5 
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The slopes of the momentum streamlines in Fig. 5-5 are 
nearly equal to the slopes of the mass flow streamlines at the 
same position for Z>7 inches, indicating negligible momentum 
is added after this axial position. Nearer the anode, however, 
the momentum streamline slopes are much less than the corres- 
ponding mass streamline slopes (i.e., the mass in a given stream 
line gains momentum) . This shows that even after the jet has 
left the anode orifice, there is an appreciable increase in 
momentum, especially near the centerline. This observation 
is in agreement with and derives from the observed velocity 
increase along the centerline (Ref. 5-1) . If this increase 
of momentum were strictly a consequence of the expansion of 
the jet into a vacuum, the momentum of the central portion 
would increase at the expense of the momentum content of the 
outer regions so the edge momentum streamlines would have a 
greater slope than the edge mass streamlines. Figures 5-3 
and 5-5 do not show this property. They show an increase in 
momentum of the entire jet of about 8N or 10%, suggesting 
small fringing magnetic fields are still accelerating the 
plasma well after it leaves the anode orifice. It must be 
recalled, however, that the momentum flux measured near the 
jet edges varied over a range of 100% (Sec. 3. A. 2), so the 
precise location and slope of the streamlines near the jet 
edges is uncertain. The precision of the data near the edges 
is not sufficient to decide between vacuum expansion or down- 
stream acceleration as an explanation of the observed center- 
line velocity increase. 

The total axial momentum (thrust) of 79 newtons used in 
Fig. 5-5 is the sum of the measured electrothermal thrust, 

(T . as reported in Chap. 4) , and the predicted electromagnetic 
thrust (T , Ref. 5-3) . It is assumed that these two forces 
dominate. The close agreement between this applied force (79 N) 
and the total integrated axial momentum in the beam (76 N) is a 
strong argument for all the assumptions used to calculate these 
numbers and also for the accuracy of the measurements. This 
argument is made even stronger by considering the predicted and 
actual dependance of thrust on current and mass flow rate. 
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Figure 5-6 shows the momentum balance obtained at the 

other mass flow rates and arc currents investigated (Table 

II). The barred verticle lines are reproduced from Fig. 3-7 

to show the measured total axial momentum content of the jet 

and are labeled with the nominal mass flow rate. The lowest 

dashed line shows the predicted T (I) which is common for all 

mass flow rates, and the other lines show the sum of this 

T (I) plus the equation developed in Chap. 4. In all except 
Gin 

one case, the total axial momentum derived from the flux pro- 
files equals the T (I) + T . ( I , m) within the accuracy of the 
data. This balance over most of the range of parameters in- 
vestigated, gives added confidence in the techniques used to 
obtain these numbers (see Chapters 3 and 4). At the lower 
currents the sensitivity limit of the pressure probe was 
approached, and back calculations showed an easily explained 
(Chap. 2) systematic error of one torr would reconcile any 
disagreement . 

5 . 3 Energy flow 

Figure 5-7 shows the flow pattern of the total energy of 

motion of the plasma produced by a mass flow rate of 5.9 g/sec 

and an arc current of 15.6 kA. This pattern was calculated by 

multiplying the total momentum flux (/^o u| |u| taken from 

the smoothed curves of Chap. 3) times 1/2 the component of 

the velocity normal to the area of integration (Fig. 5-2), 

times the weighting factor 2 7TR, integrating from R=0 to R and 

expressing this integral as a fraction of the total electrical 

6 

energy supplied to the arc (VI = 2.56 X 10 watts, see Fig. 

5-8) . This fraction was plotted against R, at the R positions 
corresponding to 10% through 35% cross plotted on Fig. 5-7 for 
Z = 5 and 11 inches, and the points connected to form flow 
energy streamlines. The most significant information gained by 
this process was that each element of the jet gained appreciable 
energy from Z=5 to Z=ll inches and that as much as 37% of the 
input power eventually appeared as directed kinetic energy. 
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KINETIC ENERGY STREAMLINES 

FIGURE 5-7 


In contrast to the observed mass deficit, it is no 
surprise that all of the energy does not appear as kinetic 
energy. Losses are expected to the electrode surfaces and 
in ionizing, excitation and heating the plasma. Sufficient 
hard experimental data are not yet available to pinpoint the 
location of the remaining 63% of VI = 2.56 MW of input power, 
but with what is available and a few reasonable assumptions a 
fairly accurate estimate is possible. That is: 

1) 11% or 0.31 MW are lost to the anode through the 
anode sheath voltage drop (Ref. 5-5) ; 

2) 16% or 0.43 MW or % of the cathode sheath power 
is lost to the cathode (Ref. 5-6) ; and 

3) 12% or 0.33 MW are lost in heating, exciting and 
ionizing the plasma to a l^eV equilibrium state 
(Refs. 5-6, 7,8). 
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The sum of all these estimates (40% plus' the 37% of kinetic 
energy is only 77% of the measured . input power, but this 
absence of an exact power balance is to be expected from the 
large uncertainties associated with the losses, especially 
the frozen flow losses of item 3, and losses (e.g. radiation) 
not included in the power balance. 

5.4 Performance Parameters 

Using the mass, momentum and flow energy profiles found 
above, many properties of the jet at the standard operating 
conditions may be calculated and used to discuss the possible 

mechanisms of the accelerator. The density, for instance, 

0 0 —6 3 

is fiu/u , typically 4 x 10 kg/m . Also, with reasonable 

assumptions, such as assuming most of the ion thermal energy 
has been converted to flow energy at Z = 11 inches, the prop- 
erties of the jet as it leaves the anode plane can be cal- 

4 

culated to be: ion temperature = 10 eV; Re=10 ; static pressure 
= 0.021 atm (nearly equal to the measured chamber pressure 
of 0.017 to 0.021 atm)? and, velocity, temperature and static 
pressure are nearly uniform across the jet. This in turn, 
implies the profile losses are low so that most (90%) of the 
kinetic energy (0.95 MW) contributes to useful thrust. 

Most of such calculations, however, will not be partic- 
ularly informative until the accuracy of the results is con- 
firmed by validating the velocity profile and inlet mass flow 
rate data. Even without this validation, however, it is useful 
to consider the performance of the accelerator as a space 
thruster, since the results suggest a promising direction for 
thruster development which may produce significant improve- 
ments in MPD accelerator performance. 

Since the discovery of the MPD accelerator, it has been 
known that performance calculations can be misleading. The 
most widely recognized error occurs at the so-called "starved" 
operating conditions. At low (or zero) inlet mass flow rates, 
an MPD accelerator still produces about the same thrust as with 
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higher mass flows (Chap. 3) . The mass accelerated is usually 
ingested or mass ablated from the arc chamber surfaces. Con- 
sequently, at these conditions, quasi-steady specific impulse 

0 m 

or, T/m, and thrust efficiency, T /2m, calculated on the basis 
of input mass flow rates do not represent the true performance 
of the accelerator as a space propulsion system. Table IV 
shows the performance parameters calculated on this basis, 
and the results for the "starved" condition I = 27 kA, and 
m = 1.9 g/sec illustrate the problem. An efficiency of 132% 
is impossible and a specific impulse of 11200 sec is improb- 
able. Indirect methods, such as measuring the electrode and 
insulator mass loss, have established that ablation occurs at 
high currents and/or low inlet mass flow rates ("starved con- 
ditions"), but these methods are not sufficiently precise to 
allow reliable calculations of performance (Ref. 5-10) . 

The techniques developed in this chapter can be used to 
give a more accurate performance evaluation for the " starved" 
operating conditions. Using the methods outlined in Sec. 5.1, 
the actual mass flow rate in the jet can be calculated. Using 
this mass flow rate in the formula T/m will give a much more 
realistic value for the specific impulse of a starved MPD 
accelerator. The result will still be too high, however, since 
as shown for the standard conditions in Sec. 5.1, the mass 
leakage around the jet is indetectable. 

A second way in which performance calculations such as 
Table IV are misleading, is concerned with thrust efficiency, 
T/m. The thrust efficiency is commonly thought of as the 
fraction of the total input power which appears as useful 
thrust. This is the ratio of the axially directed kinetic 
energy flow power to electrical power or 

yj = 1/2 l u z dA/lV 

If the velocity and momentum flux are uniform across the jet 
at a point sufficiently downstream so that the static pressure 
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TABLE IV 

Performance Parameters 
Based on INPUT Mass Flow Rates 
During Quasi-Steady Operation 


I 

A 

T 

T/A 

I sp 

T 2 /2m 

IV 

n* 

A 

kg/sec 

N 

m/ sec 

sec 

watts 

watts 

% 

7 . 75xl0 3 

5.9xl0" 3 

21*5 

3 3 /4 xlO 3 

3 . 8xl0 2 

0. 0403xl0 6 

. 98xl0 6 

4.11 

KO 

• 

00 

1.9 

22*5 

11.8 

12 

0.133 

1.2 

11.1 

11.7 

3.8 

43 

11.3 

11.5 

0.243 

2.05 

13.5 

15.6 

5.9 

79 

13.4 

13.7 

0.54 

2.56 

21 


3.9* 

79 

20.2 

20.7* 

0.8 

2.56 

33* 

26 

5.9 

220 

37.5 

38 

4.13 


17.4 

27 

1.9 

210 

110 

112 

11.5 

8.7 

132 

27*5 

23 

240 

10.4 

10.6 

1.04 


15.5 

36*5 

36 

450 


12.7 

2.86 

11.0 

26 


* based on mass flow rate in the jet and observed profiles. 
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is negligible, 

h dA = h 


the thrust power can be written as 

2 


/V 

) </> u 


y 

,) 


dA->- h 


f/ Bn : 


dA J 

“dA 


= T / 2m 


As shown in Sec. 5.3, however, the approximation of uni- 
formity is not valid. Both the velocity and momentum flux 
are functions of position at Z = 11 inches from the anode. 

Even in the anode plane where the momentum flux is nearly 
uniform across the jet, there is a leakage of plasma around 
the jet which does not participate in the MPD acceleration 

2 o 

process. The thrust efficiency as measured by T /2m, there- 
fore, is a measure not only of energy efficiency, but also 
gas or plasma leakage from the thruster. This is clearly 
misleading in a thruster development program since the pro- 
cedures used to improve the efficiency with which energy is 
converted into useful thrust are clearly different from the 
procedures required to stop a plasma or gas leak. Using the 
measured momentum flux and velocity profiles in the 'integral 
relation above, the efficiency with which energy is converted 
to thrust for the standard conditions was calculated to be 
33%. This is over half again as big as the value (21%) of 
T /2m for these conditions. Presumably, the T^/2m calculation 
is equally misleading for all the other operating conditions 
listed in Table IV as well as the reported results in Refs. 

5-4 and 5-8. The magnitude of the difference discovered for 
the standard conditions strongly suggest that the techniques 
developed in this chapter be used to separate the effect of 
leakage from MPD accelerator inefficiencies in any thruster 
optimization program. 

5.5 Thruster Development Program 

Finally, the new information gained by profile measure- 
ments shows how performance can be improved. Previous thruster 
optimization studies at AVCO (Ref. 5—4) have shown that thrus- 
ter performance is related to the parameter I 2 /m. (The results 
of Chap. 4 give a hint as to why this should be so) . Their 
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results indicated that for tjie "starved" condition, I /m>40 
2 

kA -sec/g, erosion was a problem, and they reasoned that for 
2 # 

I /m<40, the "over-fed" condition, some of the input mass 
would bypass the acceleration region. The "standard" opera- 
ting conditions in the present experiments were originally 
chosen to be "matched" (I /m' w 40) since the AVCO results im- 
plied this operating condition optimized thruster performance 
with respect to leakage and erosion. The profile studies in 
Sec. 5-1 show that even at this "matched" condition, about 
50% of the inlet mass bypasses the acceleration region. Mass 
utilization is poor, reducing the specific impulse from a 
potential 2100 sec to 1400 sec and the thrust efficiency from 
a potential 33% to 21%. In this case, the profile measurements 
provided a more direct and accurate indication of "starved" and 
"overfed" conditions than provided by the I /m correlation and 
therefore a better basis for thruster optimization. This re- 
sult suggests two complimentary programs for improving the 
thruster performance of an MPD accelerator. These sure: 

1) Stop the mass loss through the annular space between 
the jet and the anode lip. Two methods of accom- 
plishing this suggest themselves. The anode Orifice 
size could be reduced, as it is in the AVCO thrusters 
(Ref. 5-4) . There is some indication that this geom- 
etry change may decrease the efficiency by increasing 
the terminal voltage, so it will be desirable to measure 
the velocity and momentum flux profiles for a series of 
anode orifice sizes to find the optimum geometry. Al- 
ternatively, the annular space may be blocked by an in- 
sulator. The leaking plasma velocity is apparently low 
velocity so erosion of the insulator may not be serious, 
and since the leaking plasma apparently does not par- 
ticipate in the MPD acceleration process, the insulator 
may not degrade efficiency. The key point in using 
this or similar techniques for improving mass utiliza- 
tion is recognizing the problem and measuring the re- 
sult of any geometry modification. 
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2) The goal of the second program would be to find the 
optimum operating conditions for a given geometry. 

This would require measurement of the velocity and 
momentum flux profiles (as well ss m, I, V and erosion) 
over a range of m and I. Using the techniques out- 
lined in this chapter, the true specific impulse, ef- 
ficiency and mass utilization would be calculated at 
each operating condition. From this information, the 
optimum thruster operating conditions for any specified 
mission could be calculated,, 

In view of the large inefficiencies detected by the in- 
vestigation described in the preceding chapters, either or 
both of these thruster development programs could reasonably 
be expected to produce a major improvement in MPD thruster 
performance. 
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